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We present a theoretical survey of magnetocrystalline anisotropies in (Ga,Mn)As epilayers and 
compare the calculations to available experimental data. Our model is based on an envelope function 
description of the valence band holes and a spin representation for their kinetic-exchange interaction 
with localised electrons on Mn^"^ ions, treated in the mean-field approximation. For epilayers with 
growth induced lattice-matching strains we study in-plane to out-of-plane easy-axis reorientations as 
a function of Mn local-moment concentration, hole concentration, and temperature. Next we focus 
on the competition of in-plane cubic and uniaxial anisotropies. We add an in-plane shear strain 
to the effective Hamiltonian in order to capture measured data in bare, unpatterned epilayers, and 
we provide microscopic justification for this approach. The model is then extended by an in-plane 
uniaxial strain and used to directly describe experiments with strains controlled by postgrowth 
lithography or attaching a piezo stressor. The calculated easy-axis directions and anisotropy fields 
are in semiquantitative agreement with experiment in a wide parameter range. 



I. INTRODUCTION 

Dilute moment ferromagnetic semiconductors, such as 
(Ga,Mn)As, are particularly favourable systems for the 
research in basic spintronics phenomena and towards 
potential applications in memory and information pro- 
cessing technologies. For typical doping levels 1-10% of 
Mn the magnetic dipole interactions and corresponding 
shape anisotropies are 10-100 times weaker in (Ga,Mn)As 
than in conventional dense-moment ferromagnets. Con- 
sequently, magnetocrystalline anisotropy plays a decisive 
role in the process of magnetisation reversal. Despite the 
low saturation magnetisation the magnetic anisotropy 
fields reach ~ 10-lOOmT due to the large spin-orbit cou- 
pling. 

The dependence of magnetic properties of (Ga,Mn)As 
epilayers on doping, external electric fields, temperature, 
and on strain has been explained by means of an effective 
model of Mn local moments anti-ferromagnetically cou- 
pled to valence band hole spins. The virtual crystal k • p 
approximation for hole states and mean-field treatment 
of their exchange interaction with Mn d-shell moments 
allow for efficient numerical simulationsJ^'^'^ The ap- 
proach has proved useful in researching many thermody- 
namic and magneto-transport properties of (Ga,Mn)As 
samples with metallic conductivities,— such as the mea- 
sured transition temperatures^^i^ii!^ the anomalous Hall 
effect (^iiSiiiiii^ anlsotroplc magneto resistancej^iiiii^ii^iiii 
spin-stiffness,— ferromagnetic domain wall widths,— 
Gilbert damping coefficient )^^i-^^ and magneto-optical 
coefficients In this study we systematically ex- 
plore the reliability of the effective model in predicting 
the magnetocrystalline anisotropies of (Ga,Mn)As epi- 
layer and micro-devices. In our comparisons to experi- 
ment we include an extensive collection of available pub- 
lished and unpublished measured data. 

Sec. |lT] reviews key elements of the physical model of 
(Ga,Mn)As and of the corresponding effective Hamilto- 



nian used in our study. Special attention is given to 
mechanisms breaking the cubic symmetry of an ideal 
zinc-blende (Ga,Mn)As crystal. The lattice mismatch be- 
tween the epilayer and the substrate, producing a growth- 
direction strain, is responsible for the broken symmetry 
between in-plane and out-of-plane cubic axes. Micro- 
scopic mechanism which breaks the remaining in-plane 
square symmetry in unpatterned epilayers is not fully 
understood. However, it can be modelled by introducing 
an additional uniaxial in-plane strain in the Hamiltonian. 
In Sec. Ill Al we discuss the correspondence of this effec- 
tive approach and a generic k • p Hamiltonian with the 
lowered symmetry of the p-orbital states which form the 
top of the spin-orbit coupled valence band. Sec. lIIBl pro- 
vides brief estimates of the shape anisotropy in thin-film 
(Ga,Mn)As epilayers and micro(nano)-bar devices. 

Sections Uni and |IV] give the survey and analysis of 
theoretical and experimental data over a wide range of 
strains, Mn moment concentrations, hole densities, and 
temperatures. Sec. IIIlAI focuses on the easy-axis switch- 
ing between the in-plane and out-of-plane directions. 
Sec. nil Bl studies the competition of cubic and uniax- 
ial in-plane anisotropies. Sec. IIII CI provides compari- 
son based on anisotropy fields extracted by fitting the 
calculated and experimental data to the phenomenologi- 
cal formula for the magnetic anisotropy energy. Sec. IIVI 
studies in-plane easy axis reorientations in systems with 
additional in-plane uniaxial strain introduced experimen- 
tally by post-growth treatment of epilayers. Finally, in 
Sec. |V] we draw conclusions and discuss the limitations 
of our theoretical understanding of magnetic anisotropies 
in (Ga,Mn)As. 



II. MAGNETIC ANISOTROPY MODELLING 

We use the effective Hamiltonian approach to calculate 
the magneto-crystalline anisotropy energy of a system of 
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itinerant carriers exchange coupled to Mn local moments. 
The k • p approximation is well suited for the descrip- 
tion of hole states near the top of the valence band in a 
(III,Mn)V semiconductor. The strong spin-orbit interac- 
tion makes the band structure sensitive to the direction 
of the magnetisation. The Hamiltonian reads: 



in more detail. Small deformation of the crystal lattice 
can be described by a transformation of coordinates: 



n = Hkl + Jpd J2 • - i^i) + 



(1) 



Ti.KL is the six-band Kohn-Luttinger Hamiltoniansa in- 
cluding the spin-orbit coupling (see Appendix |^ . We 
use GaAs values for the Luttinger parameters. Hstr 
is the strain Hamiltonian discussed in the following sec- 
tion. The second term in Eq. ([T]) is the short-range 
antiferromagnetic kinetic-exchange interaction between 
localised spin S/ {S ~ 5/2) on the Mn^+ ions and 
the itinerant hole spin s, parametrised by a constant^'^ 
Jpd = 55 meVm~^. In the mean- field approximation it 
becomes Jpd^Mn {S) M • s. The explicit form of the 6x6 
spin matrices s is given in Ref. 0. M is the magneti- 
sation unit vector and Nmu = 4a;/aQ is the concentra- 
tion of Mn atoms in Gai-xMn^jAs (ao is the lattice con- 
stant). Note that the Fermi temperature in the studied 
systems is much higher than the Curie temperature so the 
smearing of Fermi-Dirac distribution function is negligi- 
ble. Therefore, finite temperature enters our model only 
in the form of decreasing the magnitude of magnetisa- 
tion |M| = SBs{Jpd{s)/kBT), where Bs is the Brillouin 
function, (s) is the hole spin-density calculated from the 
mean- field form of Eq. ^ . 

We emphasise that the above model description is 
based on the canonical Schrieffer-Wolf transformation of 
the many-body Anderson Hamiltonian. For (Ga,Mn)As 
the transformation replaces the microscopic hybridisa- 
tion of Mn d-orbitals with As and Ga sp-orbitals by the 
effective spin-spin kinetic-exchange interaction of L = 
0,5 = 5/2 local Mn- moments with host valence band 
statesi^ Therefore, the local moments in the effective 
model carry zero spin-orbit interaction and the magneto- 
crystalline anisotropy is entirely due to the spin-orbit 
coupled valence-band holes. The M-dependent total en- 
ergy density, which determines the magneto-crystalline 
anisotropy, is calculated by summing one-particle ener- 
gies for all occupied hole states in the valence band, 

Etot{M) Sn(k,M)/(i;„(k,M))d3A;, (2) 

n=l 

where 1 < m < 6 is the number of occupied bands 
/(£'„(k)) is the Fermi distribution function at zero tem- 
perature. 



A. Beyond the cubic symmetry of the GaAs host 

The k • p method provides straightforward means of 
incorporating elastic strains j^ i^"^'^^ which we now discuss 







(3) 



where is the strain tensor. Expressing Hkl hi 
r' coordinates leads to extra terms dependent on the 
strain that can be treated perturbatively. The resulting 
strain Hamiltonian has the same structure as the Kohn- 
Luttinger Hamiltonian with kikj replaced by . (For de- 



tailed description of Hstr see Eq. (|B2p in the Appendix.) 

Lattice matching strain induced by the epitaxial 
growth breaks the symmetry between in-plane and out- 
of-plane cubic axes. Corresponding non-zero components 
of the strain tensor read Cxx = Syy = £o = ~'§^^zz = 
[as — ao)/ao where as and ao are the lattice constant of 
the substrate and the relaxed epilayer, respectively, and 
Ci2, cii are the elastic moduli. Typical magnitudes are 
eo - 10-4 - 10-2. 

As we discuss in Sec. IIVI relaxing the growth strain 
in microbars in transverse direction produces a uniaxial 
symmetry breaking in the plane, described by a combina- 
tion of exx 7^ Syy and e^y strains, depending on the crys- 
tal orientation of the microbarj ^i^^i^^i^^ The magnitudes 
range between zero and the growth strain. Additional in- 
plane uniaxial anisotropy effects can be also induced by 
piezo stressors. The typical magnitude achieved 
by commercial stressors^ at low temperature is of the 
order of 10~*. 

An unpatterned bulk (Ga,Mn)As epilayer can also 
show broken in-plane symmetry, most frequently between 



th e [110] and 110] directions (see e.g. Refs. 

mililiolillliSliSi])- For convenience and for direct com- 
parison with effects mentioned in the previous paragraph 
we model this "intrinsic" in-plane uniaxial anisotropy by 
e™*. We fix its sign and magnitude for a given wafer by 
fitting to the corresponding measured anisotropy coeffi- 
cients. To narrow down the number of fitted values for 
e™* in the extensive set of experimental data which we 
analyse, we assume that e™* describes effectively a sym- 
metry breaking mechanism induced during growth and 
its value does not change upon the post-growth treat- 
ments, including annealing, hydrogenation, lithography 
or piezo-stressing. 

We point out that an in-plane strain has not 
been detected experimentally in the bare unpatterned 
(Ga,Mn)As epilayers. It is indeed unlikely to occur as 
the substrate imposes the cubic symmetry. The possi- 
bility of transfer of the shear strain from the substrate 
to the epilayer was ruled out by the following test ex- 
periment. A 50 nni (Ga,Mn)As film was grown on GaAs 
substrate. An identical film was grown on the opposite 
side of the neighbouring part of the same substrate. Both 
samples developed uniaxial magnetic anisotropy along a 
diagonal but the easy axes were orthogonal to each other. 
If there were a uniaxial strain in the substrate responsi- 
ble for the uniaxial anisotropy in the epilayer, the easy 
axes in the two samples would be collinear. Nevertheless, 
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wliere tlie diagonal terms of the unperturbed 6-band 
Hamiltonian corresponding to atomic orbital levels are 
set to zero. The symmetries of the tetrahedron (zinc- 
blend) point group Td narrow down the number of non- 
vanishing independent matrix elements, represented by 
Kohn-Luttinger parameters. The summation over neigh- 
bouring energy levels runs only through the Fi and r4 
states of the conduction band as other levels are excluded 
due to the parity of the wave functions or by the large 
separation in energy. After including the spin-orbit in- 
teraction and transforming to a basis of total momen- 
tum eigen-states we obtain the Hamiltonian Hkl (see 
Eqs. (|A9|) and (jAlOp in the Appendix) with three in- 
dependent Luttinger parameters 71, 72, and 73, plus a 
spin-orbit splitting parameters Ag 
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FIG. 1; (Color online) Modification of the originally cubic 
in-plane magnetic anisotropy by adding a uniaxial anisotropy 
due to the shear strain e^y or due to the local potential V — 
xy£_. eo = -0.3%, p = 3 x 10^° cm-^ x = 3%, 74 is the 
additional Luttinger parameter resulting from the in-plane 
symmetry lowering and 72 is one of the Luttinger parameters 
for GaAs (see text and Eq. (|A12|| in the Appendix). 



we argue below that the effective modelling via e™* pro- 
vides a meaningful description of the "intrinsic" uniaxial 
anisotropy. 

We compare the effective Hamiltonian corresponding 
to the e™* strain with a k • p Hamiltonian in which, 
without introducing the macroscopic lattice distortion, 
the [TlO]/[110] symmetry is broken. In the derivation of 
the 6-band Kohn-Luttinger Hamiltonian originating from 
the Asp-orbitals (denoted by \X), \Y), and \Z)), the k-p 
term is treated perturbatively to second order: 



■ E 

lf{X.Y,Z} 



{t\k-p\l){l\k-p\j) 

Ei — El 



(4) 



If the tetrahedral symmetry of the GaAs lattice is bro- 
ken the number of independent parameters increases. Let 
us consider a perturbation to the crystal potential that 
removes two of the C'2 elements of group Td (rotations by 
180° about the [100] and [010] axes). The corresponding 
potential takes a form V — xy^, which mixes the Fi 
and F4(z) states of the conduction band considered in 
the summation in Eq. ^ and leaves F4(a;) and F4(?;) 
states unchanged. (^ is a fast decreasing radial func- 
tion.) Such inter- mixing of surrounding states represents 
the local symmetry lowering of the environment of the 
valence band p-orbitals. The summation over the per- 
turbed states, aTi+l3T4{z), — f3T i + aT ^{z) , F4(x), ^4:{y) 
in Eq.|4] gives rise to extra terms in the Hamiltonian Tikp- 
(The original form Tikp is given in Eq. (|A2[) in the spin 
degenerate basis listed by Eq. (|A1|) in the Appendix.) 
Assuming a weak local potential V , a » (3, we can ne- 
glect terms of quadratic and higher order dependence on 
V and obtain: 
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+ D{kl 
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Akl 



Ck.ky + D{kl + fc2) 

B{kl+kl) + 2Dk^ky 

Ckzki, 



Ckx kz 
Cky kz 
Akl + B{kl 



(5) 



where 

D^{X\py\T^{z)){TMX). (6) 

See Eq. (jASp in the Appendix giving the full expression 
for the paramater D. Elements containing the parameter 
D change the dependence of the original Kohn-Luttinger 
Hamiltonian on the fc-vector. After considering the spin- 
orbit coupling we find that the original Kohn-Luttinger 
Hamiltonian with Tistr corresponding to e™* has the 
same form as the corrected Kohn-Luttinger Hamiltonian 
TiKL with the microscopic symmetry breaking potential 



V included if we neglect the contribution of this potential 
to the diagonal elements and replace the term D(fc^ -I- ky) 
by a constant term proportional to Cxy 

Fig. [1] illustrates that the in-plane anisotropy energy 
profile due to the local potential V can indeed be ac- 
curately obtained by the mapping on the effective shear 
strain Hamiltonian. For the particular set of material 
parameters and e™* = 0.01% considered in Fig. [Tl the 
new Luttinger parameter 74 « 72/IOO, where 74 — 
-2DmQ/3h'^ (see Eq. (|A12I) in the Appendix for the def- 
inition of 72 and the other Luttinger parameters). As we 
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discuss in the following section, effective modelling us- 
ing the strain Hamiltonian with the constant el"* term is 
sufficient to capture semiquantitatively many of the ob- 
served experimental trends. Here we have demonstarted, 
that the model effectively describes a microscopic sym- 
metry breaking mechanism yielding quantitatively the 
same in-plane anisotropy energy profiles without the pre- 
sumption of a macroscopic lattice distortion. 

B. Shape anisotropy evaluation 

We conclude this theoretical modelling section by 
briefly discussing the role of shape anisotropy in 
(Ga,Mn)As thin films and microstructures. Magnetic 
shape anisotropy is due to the long range dipolar interac- 
tion. Surface divergence of magnetisation M gives rise to 
demagnetising field H^(M, r). In homogeneously mag- 
netised bodies of general shape the demagnetising field 
is a function of magnetisation magnitude and direction 
with respect to the sample. In ellipsoidal bodies the func- 
tion becomes linear in M and H^(M) is uniform in the 
body: 

iJf(M) = -^iV,,M,. (7) 



Tensor Nij is the so called demagnetising factor. In rect- 
angular prisms the linear formula ([7]) is a good approxi- 
mation and the non- uniform demagnetising factor can be 
replaced by its spatial average. For the magnetostatic en- 
ergy density of a homogeneously magnetised rectangular 
prism we get: 

E''{m) = -]^^iY,N^J{a,h,c)M,M,, (8) 

where we assume a prism extending over the volume 
—a < X < a, —b < y < b and — c < z < c in a Carte- 
sian coordinate system. Ref. [1^ shows the expression 
for Nij (a, b, c) in such prism. 

Fig. [5] shows the calculated shape anisotropy energy 
Ea = E^{Mi) - E^{M2) for a (i) thin film with a = 6 > 
c and with magnetisation out-of-plane or in-plane (Mi = 
(0,0, M), M2 = (M,0,0)), and (ii) for a bar with a > 
b ^ c and with magnetisation in-plane (Mi = (0, M, 0), 
M2 = (M, 0, 0)). In the former case the shape anisotropy 
favours in-plane easy-axis direction while in the latter 
case the easy-axis tends to align along the bar. 

As a result of the relatively low saturation magnetisa- 
tion of the dilute magnetic semiconductor, the in-plane 
vs. out-of-plane shape anisotropy Ea is only about 1.4 
kJ/m^ (0.06 T) for Mn doping x = 5% and c < a/100. 
This is in agreement with the limit of infinite 2D sheet, 
where the formula for shape anisotropy energy per unit 
volume simplifies to Ea = cos^ 9. 9 is the an- 

gle that the saturation magnetisation M subtends to the 
plane normal. The in-plane anisotropy of a bar is even 
weaker and decreases with relative widening of the bar. 




Ratio r 

FIG. 2; Shape anisotropy Ea = £"(Mi) - E'^ {Kh) of a film 
of a thickness c and a long bar of length a and width 6 as a 
function of the dimension-less ratio r as defined in the caption. 
The curves were obtained using the demagnetising factor ap- 
proximation of Ref. [45! for |M| = 0.06T which corresponds 
to Mn doping of a; = 5% at T = OK. 



In general, the shape anisotropies in the (Ga,Mn)As 
dilute-moment ferromagnet are weak compared to 
the spin-orbit coupling induced magneto-crystalline 
anisotropies and can be often neglected. 



III. MAGNETIC EASY AXES IN 
UNPATTERNED SAMPLES 

A large amount of experimental data on magnetic 
anisotropy in (Ga,Mn)As has accumulated over the past 
years. Comparison of these results with predictions of 
the effective Hamiltonian model is not straightforward 
due to the presence of unintentional compensating de- 
fects in (Ga,Mn)As epilayers. Most importantly, a frac- 
tion of Mn is incorporated in interstitial positions. These 
impurities tend to form pairs with Mnca acceptors in as- 
grown systems with approximately zero net moment of 
the pair, resulting in an effective local-moment doping 
Xeff = Xs — Xi.^ Here Xs and Xi are partial concentra- 
tions of substitutional and interstitial Mn, respectively. 
In as-grown materials, the partial concentration Xi in- 
creases with the total Mn concentration, xtot = Xs + Xi. 
For Xtot > 1.5%, dxi/dx « 0.2.^ We emphasise that in 
theory the Mn local moment doping labelled as " x" corre- 
sponds to the density of uncompensated local moments, 
i.e., to Xeff in the notation used above. Mn doping "x" 
quoted in experimental works refers typically to the total 
nominal Mn doping, i.e., to Xtot- When comparing theory 
and experiment this distinction has to be considered. 

Although interstitial Mn can be removed by low- 
temperature annealing, Xeff will remain smaller than 
the total nominal Mn doping. The interstitial Mn im- 
purities are double donors. Assuming no other sources 
of charge compensation the hole density is given by 
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The concentration of ferromagnetically ordered Mn lo- 
cal moments and holes is not accurately controlled dur- 
ing growth or determined post growth/ We acknowledge 
this uncertainty when comparing available magnetome- 
try results with theory. Throughout the paper we test 
the relevance of our model over a wide parameter range, 
focusing on general trends rather than on matching re- 
sults directly based on the material parameters assumed 
in the experimental papers. 



A. In-plane vs. out-of-plane magnetic easy axis 

In this section we study the switching between in-plane 
and perpendicular-to-plane directions of the magnetic 
easy axis. (Anisotropics within the growth plain of a 
sample are studied in Sec. IIII Bl ) Early experiments were 
suggesting that the in-plane vs. perpendicular-to-plane 
easy axis direction is determined exclusively by the sign 
of the growth induced strain in the sample. The in-plane 
easy axis (lEA) develops for compressive growth strain 
Co — (fls — ao)/ao < 0. Tensile growth strain, eo > 0, re- 
sults in the perpendicular-to-plane easy axis (PEA) . This 
simple picture was subsequently correcte d by expe rimen- 
tal results reported for example in Refs. [33k^l47E8ll49| | . 
Sign changes in the magnetic anisotropy for the same 
sign of the growth strain were observed with varying Mn 
concentration, hole density, and temperature. 

An overview of theoretical easy axis reorientations 
driven by changes of the material parameters is given 
in Figs.[3]-[6l In the plots we show the difference AE be- 
tween total hole energy density for the magnetisation ly- 
ing in-plane {Etot{M\\)) and out of plane {Etot{M±)) as a 
function of the hole density and temperature. {Etot{M\\) 
is always the smaller of Etot for magnetisation along the 
[100] and the [110] axis.) We include calculations for four 
Mn local moment concentrations to facilitate the compar- 
ison with experimental data of different nominal Mn con- 
centrations and different degree of annealing, which also 
increases the number of uncompensated local moments as 
discussed above. We note that the calculated magneto- 
crystalline anisotropics are almost precisely linear in the 
growth strain and therefore the boundaries between lEA 
and PEA in the Figs. [3]- [5] depend only very weakly on the 
magnitude of the growth strain, certainly up to the typ- 
ical experimental values |eo| < 1%. Magneto-crystalline 
anisotropy diagrams presented in this section for a com- 
pressive strain Cq = —0.2% are therefore generic for all 
typical strains, with the lEA and PEA switching places 
for tensile strain. 

Solid arrows in Figs. [3] -[6] mark easy-axis behaviour as 
a function of temperature and doping that has been ob- 
served experimentally. The dashed arrows correspond to 
theoretical anisotropy variations that have not been ob- 
served experimentally. At low hole densities, increasing 
temperature (marked by arrow (1)) induces a reorienta- 
tion of the easy axis from a perpendicular-to-plane to 
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FIG. 3: (Color online) Anisotropy energy AE = E{M\\ 
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calculated for x — 



eo = -0.2%, 



0. Positive(negative) AE corresponds to lEA(PEA). Arrows 
mark anisotropy transitions driven by change of temperature 
or hole density. 
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FIG. 4: (Color online) Anisotropy energy AE = E{M\i) - 
E{M±) [kJm-^] calculated for x = 6%, eo = -0.2%, e^y = 
0. Positive(negative) AE corresponds to lEA(PEA). Arrows 
mark anisotropy transitions driven by change of temperature 
or hole density. 



an in-plane direction. With decreasing x this transition 
shifts to lower hole densities; at a; = 2% the theoretical 
densities allowing for such a transition reach unrealisti- 
cally low values for a ferromagnetic (Ga,Mn)As mate- 
rial with metallic conduction. Warming up the partially 
compensated samples (marked by arrow (2)) has no re- 
orientation effect and the easy axis stays in-plane. There 
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Temperature [T/TJ ° Temperature [T/TJ 



FIG. 5: (Color online) Anisotropy energy AE — E(M^\) — 
E{M_l) [kJm"^] calculated for x = 4%, eo = -0.2%, e^y = 
0. Positive (negative) AE corresponds to lEA(PEA). Arrows 
mark anisotropy transitions driven by change of temperature 
or hole density. 



FIG. 6: (Color online) Anisotropy energy AE — E{M\^) — 
E{Mi_) [kJm"^] calculated for x = 2%, eo = -0.2%, e^y = 
0. Positive(negative) AE corresponds to lEA(PEA). Arrows 
mark anisotropy transitions driven by change of temperature 
or hole density. 



are no exceptions to this behaviour at different Mn con- 
centrations. Finally, increasing temperature of a very 
weakly compensated (fully annealed) sample can cause 
switching of the theoretical easy direction from in-plane 
to perpendicular-to-plane (marked by arrow(3)), with the 
exception of the low Mn concentrations. 

The techniques used to increase the hole density in 
the experimental works discussed in this section are the 
postgrowth sample annealing and annealing followed by 
hydrogen passivation/depassivation.'*^ The latter method 
yields solely a change of hole density, whereas the former 
is associated also with an increase of the effective Mn 
concentration and a decrease of the growth strain. The 
growth strain is caused to a large extent by Mn atoms in 
interstitial positions,™ which are removed by the anneal- 
ing. The simultaneous increase of hole density and effec- 
tive Mn concentration due to annealing implies a transfer 
between the phase diagrams of Figs. [3] -[S] accompanying 
the transitions marked by arrows (4) - (6). We argue that 
the remarkable similarity of the four diagrams assures a 
meaningfull qualitative comparison with the effect of an- 
nealing even within a given diagram. 

We now discuss individual measurements and compare 
with theoretical diagrams in Figs. [3]- [S] Ref. ^] reports 
experiments in a 50 nm thick (Ga,Mn)As epilayer nomi- 
nally doped to a; = 6—7% and grown on a GaAs substrate 
under compressive strain. The sample is first annealed to 
lower the number of interstitial Mn, then hydrogenated 
to passivate virtually all itinerant holes and finally de- 
passivated in subsequent steps by annealing. The hole 
density was not measured but for the given Mn doping 
we expect the density in the range of p ~ 10^" — lO^^cm"'^ 



after depassivation. The low temperature (T = 4 K) re- 
orientation from PEA to lEA induced by successive de- 
passivations and detected indirectly by anomalous Hall 
effect measurement in Ref. [i^ matches the transition 
marked by arrow (4) in Figs. [3] -[H 

Magnetic h yste resis loops measured by the Hall resis- 
tivity in Ref. [49| reveal easy axis reorientations induced 
by annealing or increasing temperature in material with 
nominal Mn doping x = 7%. This (Ga,Mn)As epilayer 
was grown on a (In,Ga)As buffer which leads to a ten- 
sile strain. (Recall that the anisotropy energy AE is an 
odd function of the growth strain so the lEA and PEA 
regions have to be interchanged in Figs. OlS] when con- 
sidering tensile strain.) Again, the hole density is not 
known and can be estimated to p ~ 10^° — lO^^cm"'^. 
After annealing, the material exhibits perpendicular-to- 
plane easy axis at 4 K and no reorientation occurs during 
heating up to 115 K (Tc « 120-130 K in this material). 
Such behaviour corresponds to arrow (2) of Fig. [3]or Fig. 
|3l The as-grown sample has lEA at 4 K and PEA at 
22 K. This easy axis reorientation corresponds to arrow 
(1), again considering a tensile strain. The as-grown and 
annealed samples both share PEA at elevated temper- 
ature. Such a stability of the easy axis while changing 
the hole density corresponds to arrow (5). Theoretical 
anisotropy variations described by arrows (3) and (6) are 
not observed in Ref. 49] 

Ref. [3^ presents measurements in compressively 
strained (Ga,Mn)As epilayers grown on a GaAs sub- 
strate. The reported nominal Mn concentrations are 
X = 5.3% and x = 3% with compressive growth strain 
eo = —0.27% and eo = —0.16%, respectively, as inferred 
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from x-ray diffraction measurement of the lattice param- 
eter. The higher doped material was partially annealed 
for several different annealing times. The hole density 
was not measured but likely increases substantially with 
annealing. The as- grown x = 5.3% sample at 5 K ex- 
hibits PEA, which changes to lEA upon warming up to 
22 K. This anisotropy variation is not observed for sam- 
ples subject to long annealing times. Such a result is 
consistent with Ref. ^ and corresponds to the theoret- 
ical predictions marked by arrows (1) and (2) of Fig. [S] 
for increasing temperature of the as-grown and annealed 
sample, respectively. Again, the effect of annealing is in 
good agreement with anisotropy behaviour predicted for 
low (high) temperature represented by arrow (4) (arrow 
(5)), however, there is no experimental counterpart of 
transitions marked by arrows (3) and (6). The sample 
doped to a; = 3% was not annealed and no transition 
from PEA to lEA is observed upon warming. The be- 
haviour corresponds to arrow (2) in Fig. [S]or[5l 

Ref. [Slj already reports a successful compari- 
son of measured magnetic anisotropy and theoretical 
predictionSfi Among other samples, it presents a com- 
pressively strained (Ga,Mn)As epilayer with nominal Mn 
concentration x = 2.3% (inferred from x-ray diffraction 
measurement). A superconducting quantum interference 
device (SQUID) measurement of this as-grown sample 
shows PEA at 5 K and lEA at 25 K, corresponding to 
anisotropy variation marked by arrow (1) in Fig. [5] (oc- 
curring only for a very narrow hole density interval). 

Ref. [12] presents (Ga,Mn)As epilayers with compres- 
sive and tensile strain grown on GaAs and (In,Ga)As 
buffers, respectively, with nominal Mn concentration 
X = 3% inferred from reflection high energy electron 
diffraction (RHEED) oscillations measured during the 
molecular-beam epitaxy (MBE) growth. Two of the sam- 
ples are annealed and magnetic anisotropy is investigated 
at 5 K. The tensile strained sample has its easy axis 
aligned perpendicular to the growth plane and the com- 
pressively strained sample has an in-plane easy axis. This 
observation is in good agreement with our theoretical 
modelling. 

Finally, Ref. [46| shows a transition from PEA to IMA 
upon increasing temperature or change of hole concen- 
tration (induced by gating in this case) . The sample is a 
(In,Mn)As epilayer grown on an InAs, and its magnetic 
anisotropy is described consistently by our model when 
the appropriate band parameters are used. 



B. In-plane anisotropy: Competition of cubic and 
uniaxial components 

As we discussed in the previous section, the magnetic 
easy axis(axes) is in the plane of (Ga,Mn)As/GaAs films 
over a wide range of dopings. Experimental works in bare 
(Ga,Mn)As epilayers discussed in this section show that 
the in-plane magnetic anisotropy has cubic and uniax- 
ial components. Typically, the strongest uniaxial term is 
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FIG. 7: (Color online) Magnetic anisotropy energy AE = 
— -E[ioo] as a function of the in-plane magnetisation orien- 
tation M — |M| [cos (j}, sin 0, 0] and its dependence on material 
parameters. Magnetic easy axes (marked by arrows) change 
their direction upon change of hole density p given in units 
u = 10^" cm~'^ at Mn local moment concentration x — 5%, 
shear strain Cxy = 0, and zero temperature. 

along the in-plane diagonal ([110] /[1 10]) direction. (A 
weak uniaxial component along the main crystal axes 
([100]/ [010]) has also been detected 4^ii^) The theoretical 
model used so far to describe the easy axis reorientation 
between the in-plane and out-of-plane alignment, assum- 
ing the growth strain, can account only for the cubic 
in-plane anisotropy component. In this case we find two 
easy axes perpendicular to each other either along the 
main crystal axes or along the diagonals depending on 
the Mn concentration and hole density, as shown in Fig. 
[T] In order to account for the uniaxial component of the 
in-plane [110]/[110] anisotropy in bare (Ga,Mn)As epilay- 
ers the elastic shear strain Cxy is incorporated into our 
model as discussed in Sec. HT] (For brevity we omit the 
index "int" in the following text and reintroduce the in- 
dex only when additional real in-plane strains are present 
due to micro-patterning or attached piezo-stressors.) The 
superposition of the two components results in a rich phe- 
nomenology of magnetic easy axis alignments as reviewed 
in Fig.[5]-[ini 

Fig. [8] shows an example with easy axes aligned close to 
the main crystal axes [100] and [010] at Mn local moment 
concentration x — 5%, hole density p = 3 x lO^'^cm^'^, 
and a weak shear strain exy = 0.01%. For a stronger 
shear strain Cxy = 0.03% the cubic anisotropy is no longer 
dominant and the easy axes "rotate" symmetrically to- 
wards the diagonal [110] direction until they merge for 
Gxy ^ 0.05%. As explained in detail in Sec. |TT1 the mag- 
nitude and sign of the intrinsic shear strain Cxy enter as 
free parameters when modelling in-plane anisotropies of 
bare epilayers. 

The relative strength of uniaxial and cubic anisotropy 
terms depends also on the hole density and Mn concen- 
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FIG. 8: (Color online) Magnetic anisotropy energy Ai5 = 
— -E[iooi ^ ^ function of the in-plane magnetisation ori- 
entation M = |Ml[cos 0, sin 0, 0] and its dependence on ma- 
terial parameters. Magnetic easy axes (marked by arrows) 
change their direction upon change of magnitude of shear 
strain Cxy > at Mn local moment concentration x = 5%, 
hole density p — 3 x 10^" cm~'^, and zero temperature. 



FIG. 10: (Color online) Magnetic anisotropy energy AE = 
E^ — -Efioo] as a function of the in-plane magnetisation orien- 
tation M = |M| [cos 0, sin <p, 0] and its dependence on material 
parameters. Magnetic easy axes (marked by arrows) change 
their direction upon change of Mn local moment concentration 
X at hole density p = 3 x 10^" cm"'', shear strain Cxy = 0.01%, 
and zero temperature. 




FIG. 9: (Color online) Magnetic anisotropy energy Ai5 = 
Etf, — -E[ioo] as a function of the in-plane magnetisation orien- 
tation M = |M| [cos (j), sin 4>, 0] and its dependence on material 
parameters. Magnetic easy axes (marked by arrows) change 
their direction upon change of hole density p given in units 
u = lO^'' cm"'', at Mn local moment concentration x = 3%, 
shear strain e^y = 0.01%, and zero temperature. 



tration as shown by Fig. [9] and [TOl respectively. Both 
anisotropies are non-monotonous functions of x and p, 
compared to the linear dependence of uniaxial anisotropy 
on the shear strain. We do not show explicitly the effect 
of increasing temperature which in the mean-field theory 
is equivalent to decreasing the effective Mn concentration 
while keeping the hole density constant (as explained in 
SecHH). 



We begin the comparison of theory and experiment 
by analysing experimental studies of in-plane magnetic 
anisotropy in bare samples without lithographically or 
piezo-electrically induced in-plane uniaxial strain. Ex- 
perimental results are summarised in Tab.|Tl Samples are 
identified by nominal Mn concentration and hole density 
or annealing as given by the authors. Typically, the hole 
density is in the range lO^*' - lO^^cm""^. All samples are 
thin (Ga,Mn)As epilayers deposited by MBE on a GaAs 
substrate. According to our calculations, the compres- 
sive growth strain has a negligible effect on the interplay 
of cubic and uniaxial in-plane anisotropies. 

Tab. [T] shows the largest measured projection of the 
easy axis (axes) on the main crystal directions ([100], 
[010], [110], [110]) in the corresponding sample. (Note 
that unlike our theoretical calculations of the full in-plane 
anisotropy profile, most experiments listed in Tab. |T] re- 
port only projections of the magnetisation to the main 
crystal directions. Studies using anisotropic magneto- 
resistance (AMR) to map the easy axis direction pre- 
cisely are discussed in Sec. IIII Cl and HVl ) Tab. [T] includes 
a column labelled as EAq giving the largest easy axis 
projection at low temperatures (typically 4 K) and a col- 
umn labelled as EA^c corresponding to measurements at 
temperatures close to Tc- This simplified overview of the 
temperature-dependence of the in-plane anisotropies re- 
flects the nature of available experimental data. The fer- 
romagnetic resonance (FMR) spectra are typically pro- 
vided only at one high and one low temperature. More- 
over, available SQUID data reveal at most one transi- 
tion between main crystal directions corresponding to the 
largest projection of the magnetisation in the whole tem- 
perature interval. Sample No. 25 in Tab. |T] which shows 
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two transitions is the only exception to this trend. 

From Tab. [T] we infer the following general trend in 
the experimentally observed in-plane anisotropies: At 
low temperatures the in-plane anisotropy is dominated 
by its cubic component. In most cases, this leads to two 
equivalent easy axes aligned close to [100] and [010] direc- 
tions. Only in a few samples the cubic anisotropy yields 
easy-axis directions along the [110]/[1-10] diagonals at 
low temperature. The two diagonals are not equiva- 
lent, however, due to the additional uniaxial anisotropy 



component 



40.53.54.55 



At high temperatures the uniaxial 



anisotropy dominates giving rise to only one diagonal 
easy axis. Finally we note that Refs. [Hill do not iden- 
tify the correspondence between the in-plane diagonal 
easy-axis and one of the two non-equivalent crystallo- 
graphic axes [110] and [1-10] (these measurements are 
marked as (g) in Tab. This ambiguity does not affect 
the comparison with our modelling of unpatterned bare 
films since the shear strain e^y determining which of the 
two diagonals is magnetically easier is a free effective pa- 
rameter of the theory. Possibility of error in assigning the 
two non-equivalent diagonal crystallographic axes is ac- 
knowledged by the authors of Ref. [Sj] , where switching 
roles of the diagonals makes the results consistent with 
later works of the group. 

Following the strategy for presenting experimental 
data in Tab. HI we plot in Figs. [Tl] - [16] theoretical di- 
agrams indicating crystallographic axes ([100], [110] or 
[110]) with the largest projection of magnetisation as a 
function of the hole density and temperature. The com- 
parison with experimental results in Tab. [T] is facilitated 
by numbered arrows added to the diagrams, which corre- 
spond to switchings between crystallographic directions 
with the largest projection of the easy-axis, driven by in- 
creasing temperature (horizontal arrows) and hole den- 
sity (vertical arrows). 

Figs. [TT] - [14] present diagrams for different Mn con- 
centrations and for e^y = 0.01%. Anisotropy transitions 
seen in the figures are consistent with majority of the re- 
viewed experimental works, i.e., the arrows correspond to 
the experimentally observed transitions and their place- 
ment in the diagrams is reasonably close to the relevant 
experimental parameters. Figs . [TT] - [Til also demonstrate 
how the transition from the [100] to the [110] direction 
moves to higher temperatures with increasing Mn local 
moment concentration. 

Figs. [15] and [TH] address samples where the observed 
transition cannot be modelled e-xv = 0.01%. Four of 
the low doped samples in Refs. SSIIsgI 37t are modelled by 
a weaker strain, whereas one of the highly doped samples 
in Ref. [13] is modelled by a stronger strain. 

Now we discuss in detail the theoretical diagrams in 
Figs. [TT] - [TT] and compare to individual samples from 
Table [I] referred to as TjIJ-No. Fig. [TT] maps in-plane 
magnetic anisotropy at Mn local moment concentration 
X — 3% and shear strain Cxy = 0.01%. The easy axis 
reorientation of the as-grown sample T]l]-7 corresponds to 
arrow (1) in Fig.[TT] Arrow (2) in Fig. [TT] highlights the 
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TABLE I: Experimental in-plane magneto-crystalline 
anisotropies at low temperature EAit, and high temperature 
EAhT extracted from SQUID or FMR measurements: largest 
easy axis projection along [100] and [010] axes (-I-), along 
[110] axis (\), along [110] axis (/"), and along one of the 
[110]/[110] diagonals not distinguished in the experiment 
(®). Nominal Mn concentrations x reported in experimental 
studies are rounded down to percents. Hole density p [*] is 
given in units of lO^^cm"^. If the hole density is unknown 
the as-grown and annealed samples are indicated by "ag" 
and "an", respectively. Samples are ordered according 
to Mn concentration and hole density (annealed sample 
follows the as-grown counterpart when it exists). The last 
four columns label the experimental data in a way which 
facilitates direct comparison with transitions highlighted by 
arrows in the theory Figs. [TT] - 1161 Numbers in columns 
Ap, A;t, and Kht point to corresponding theory transitions 
marked by horizontal arrows, vertical arrows at low T, and 
vertical arrows at high T, respectively. The index " indicates 
correspondence of the given arrow to modelling with negative 
value of e^y. 



finite range of hole densities for which the largest pro- 
jection of the easy-axes stays along the [100] and [010] 
directions at low temperature, consistent with the be- 
haviour of the as-grown and annealed sample dJl]-9 and 
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FIG. 11: (Color online) Theoretical hole density - temper- 
ature diagrams of crystal directions with the largest projec- 



tion of the magnetic easy axis at a; = 3%, 



0.01%, 



eo = —0.2%. Arrows mark anisotropy behaviour driven by 
change of temperature or hole density explaining experimen- 
tally observed behaviour surveyed in Tab. |T] 
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FIG. 13: (Color online) Theoretical hole density - temper- 
ature diagrams of crystal directions with the largest projec- 
tion of the magnetic easy axis at a; = 7%, e^cy = 0.01%, 
eo — —0.2%. Arrows mark anisotropy behaviour driven by 
change of temperature or hole density explaining experimen- 
tally observed behaviour surveyed in Tab. [T] 
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FIG. 12: (Color online) Theoretical hole density - temper- 
ature diagrams of crystal directions with the largest projec- 
tion of the magnetic easy axis nt x — 5%, e^y = 0.01%, 
eo — —0.2%. Arrows mark anisotropy behaviour driven by 
change of temperature or hole density explaining experimen- 
tally observed behaviour surveyed in Tab. U] 




Temperature 



FIG. 14: (Color online) Theoretical hole density - temper- 
ature diagrams of crystal directions with the largest projec- 
tion of the magnetic easy axis at a; = 9%, e^y = 0.01%, 
eo — —0.2%. Arrows mark anisotropy behaviour driven by 
change of temperature or hole density explaining experimen- 
tally observed behaviour surveyed in Tab. H] 



T]I]-10. (Note that hole densities in samples Tl^9 and 
T|I]-10 were measured by the electrochemical capacitance- 
voltage profiling.) The transition from the largest easy 
axis projection along the cube edges to the [110] diag- 
onal observed in as-grown sample l]l]-8 with increasing 
temperature has no analogy in Fig. [TT] or Fig. [T^ The 
FMR measurement does not indicate switching of the 
easy axis alignment between the diagonals at any inter- 
mediate temperature. This behaviour can be explained 
only if the opposite sign of the shear strain is used to 
model the intrinsic symmetry breaking mechanism. Then 
the easy axis transition of Tjl]-8 would correspond to ar- 



row (1) in Fig.[ni 

The behaviour of as-grown samples Tll]-ll,13,15 cor- 
responds to arrow (2) in Fig. [121 The annealed sam- 
ples T1I]-12,16 exhibit the rarely experimentally observed 
domination of uniaxial anisotropy for the whole tempera- 
ture range. This behaviour is also consistently captured 
by the theory as highlighted by arrow (3) in Fig. [T^ 
Sample T]l]-14 has a dominant cubic anisotropy preferring 
[100]/ [010] magnetisation directions at low temperature 
and the easy axis aligns closer to the [110] direction at 
high temperatures. Similarly to sample T]l]-8, this transi- 
tion has no analogy in Fig.[TT]or Fig.[T21 however, it can 
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be explained assuming that the [110]/[1-10] symmetry 
breaking mechanism has opposite sign in this material 
and therefore should be modelled by a negative value of 
the effective strain e^y- Then the easy axis transition of 
T|I]-14 would correspond to arrow (2) in Fig. [121 Another 
possibility is to assume the same sign of Cxy as for the 
above samples and associate the transition in sample Tjl]- 
14 with arrow (4) in Figs. [T51 and [HI Note, however, that 
the intermediate-temperature anisotropy state with the 
largest magnetisation projection along the [110] diagonal 
seen when following the theory trend along arrow (4) has 
not been reported in the experimental study of sample 
l|l]-14. Arrows (4)-(6) in Fig. [121 correspond to measured 
anisotropy behaviour driven by increasing hole density 
in pairs of as-grown and annealed samples T|ll-ll,12, TJH 
13,14, and Tllll5,16. 

At the upper end of the investigated effective Mn con- 
centration interval the theoretical alignment of magnetic 
easy axes is mapped by Figs.[T3land[T4l Samples lH-lT to 
l]ll20 nominally doped to x = 7% were all annealed after 
growth, passivated by hydrogen plasma, and then gradu- 
ally depassivated to achieve different hole densities (mea- 
sured by high- field Hall effect). Magnetic anisotropics 
were determined by FMR. The assignment of the in-plane 
diagonal directions to the non-equivalent [110] and [110] 
crystallographic axes is not specified in this experimen- 
tal work; recall that this ambiguity is not crucial for the 
present discussion. The transition observed in these sam- 
ples from a cubic ([100]/ [010] easy directions) dominated 
anisotropy at low temperatures to a uniaxial behaviour 
at high temperatures is captured by arrows (3) and (4) 
in Figs. [Island [141 Importantly, the depassivated higher 
hole density samples T[lll9 and T[ll20 show an additional 
switching of the easy-axis from one to the other diagonal 
direction at intermediate temperatures, consistent with 
the theoretical temperature dependence along the arrow 
(4). This double transition behaviour was also detected 
in the annealed sample T]Il-25, where the temperature 
dependent magnetisation projections were measured by 
SQUID. In this experiment it is identified that the easy- 
axis first rotates towards the [110] direction at interme- 
diate temperatures and then switches to the [110] direc- 
tion at high temperatures, consistent with the behaviour 
marked by arrow (4) in Figs. [T51 and [HI 

Samples Tlll21,22 are measured only at low tempera- 
ture. Easy axis reorientation from [100] to [110] direction 
is driven by increase of hole density, which corresponds to 
arrow (6) in Fig.[T31or[T31 The hole density was measured 
by the electrochemical capacitance-voltage method. 

In-plane anisotropics of samples with x sa 2% are mod- 
elled in Fig. 1151 To obtain the cubic anisotropy domi- 
nated region at low temperatures and a transition to the 
uniaxial behaviour at high temperatures, as observed in 
samples Tl^l to TIBG, we take for this low Mn doping 
Cxy = 0.005%. (The effective strain e^y — 0.01% would 
lead to easy axis along [110] over the entire temperature 
range and for e^y = 0.001% the cubic anisotropy region 
would extend up to very high temperatures.) Arrow (1) 




Temperature 



FIG. 15: (Color online) Theoretical hole density - temper- 
ature diagrams of crystal directions with the largest projec- 
tion of the magnetic easy axis ai x — 2%, e^y = 0.005%, 
eo ~ —0.2%. Arrows mark anisotropy behaviour driven by 
change of temperature or hole density explaining experimen- 
tally observed behaviour surveyed in Tab. [11 
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FIG. 16: (Color online) Theoretical hole density - temper- 
ature diagrams of crystal directions with the largest projec- 
tion of the magnetic easy axis at a; = 7%, e^y = 0.03%, 
eo — —0.2%. Arrows mark anisotropy behaviour driven by 
change of temperature or hole density explaining experimen- 
tally observed behaviour surveyed in Tab. [H 

in Fig. [151 corresponds to easy axis switching from the 
[100] to the [ITO] direction in samples T[ll-1,2,3. Arrows 
(2) and (3) in Fig. [TH mark the behaviour of the easy 
axis driven by increasing hole density when annealing the 
sample TjUl to obtain the sample Tl^2 at low and high 
temperature, respectively. Sample Tlll4 assumes the [110] 
diagonal always harder than the [110] diagonal. A tran- 
sition from cubic to uniaxial dominated anisotropy is ob- 
served upon increasing the temperature. This behaviour 
corresponds to arrow (1) in Figs. [TSl (The hole density 
of sample TjlS, p = 3.5 x 10^°cm~'^, was determined 
by low-temperature high-field Hall effect measurements. 



however, it was not measured for samples T]l]-1,2,4.) 

Samples Tl^5 and 1^6 have their easy axis aligned 
closer to the [100]/ [010] directions at low temperatures 
and to the [110] direction at higher temperatures, simi- 
larly to sample THj-S. The SQUID measurement of mag- 
netisation projections for the whole range of temperature 
does not indicate the easy axis alignment close to the 
[110] direction at any intermediate temperature. The 
hole density of the sample T!^5, p = 1.1 x 10^"cm~'^, 
is measured by Hall effect (at room temperature) and 
its Mn concentration is inferred from X-ray diffraction 
measurement of the lattice constant. The hole density 
of the sample TJH-G is p = 4 x 10^° cm~^ (measured by 
the electrochemical capacitance-voltage method at room 
temperature) and we estimate the Mn concentration from 
the reported critical temperature, Tc = 62 K, after an- 
nealing. The described experimental behaviour does not 
correspond to predicted anisotropy transitions for rele- 
vant hole densities, Mn local moment concentrations, and 
positive shear strain. The behaviour can be explained, 
however, if the opposite sign of the shear strain is used 
to model the intrinsic symmetry breaking mechanism at 
low Mn concentration. Then the easy axis transition of 
T|I]-5,6 would correspond to arrow (1) in Fig. fTSl 

Finally we comment on the less frequent behaviour 
observed in the annealed sample Tll]-24. While its as- 
grown counterpart T]I]-23 shows the commonly seen tran- 
sition from the cubic dominated anisotropy to the uni- 
axial anisotropy with increasing temperature, marked by 
arrow (1) in Fig. [111 the annealed material has its easy 
axis aligned close to the [110] direction over the entire 
studied temperature range. Arrow (2) in Fig. [16] pro- 
vides an interpretation of this behaviour if we increase the 
magnitude of the effective shear strain. At e„j — 0.03% 
the cubic anisotropy dominated region is already strongly 
diminished and for e^y — 0.05% it vanishes completely. 
Arrows (3) and (4) then highlight within the same dia- 
gram the consistent description of the evolution of the 
experimental anisotropics, both at low and high temper- 
atures, from the as-grown low hole density sample T|l]-23 
to the annealed high hole density sample T]l]-24. 

To summarise this section, our theoretical modelling 
provides a consistent overall picture of the rich phe- 
nomenology of magneto-crystalline anisotropics in un- 
patterned (Ga,Mn)As epilayers. Our understanding is 
limited, however, to only a semiquantitative level, ow- 
ing to the approximate nature of the mean-field kinetic- 
exchange model, ambiguities in experimental material 
parameters of the studied films, and unknown micro- 
scopic origin of the in-plane uniaxial symmetry breaking 
mechanism. We remark that the effective shear strain 
we include to phenomenologically account for the exper- 
imental [110] /[1 10] uniaxial anisotropy scales with Mn 
doping {exy — 0.005a;). It brings additional confidence 
in this modelling approach as it is most likely the incor- 
poration of Mn which breaks the cubic symmetry of the 
lattice. The magnitude of the effective strain parame- 
ter falls into the range 0.005% < Cxy < 0.05% and the 
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FIG. 17: (Color online) In-plane uniaxial anisotropy as a 



function hole density at zero temperature, Cx 



0.05%, 



and eo = calculated in this work (a) and in Ref. [35| (b). 
Curves are labelled by the valence-band spin-splitting param- 
eter Bg = JpdNMnS/6 to allow for simple comparison with 
Ref. [1^. (Bg = 4.98a: in meV and in percent, respectively.) 
Dashed intervals of the horizontal axis mark regions where 
a change of temperature (inversely proportional to Bg) can 
lead to the [TlO] [110] easy axis reorientation. 



anisotropy behaviour consistent with most experimental 
works is modelled with positive sign of Cxy 



We conclude this section by a remark on numerical 
simulations of the [110] to [110] easy axis transition per- 
formed in Ref. [33] ■ The physical model employed by 
the authors of Ref. fsHj is identical to ours, neverthe- 
less, the results of the calculations do not quantitatively 
match ours, as illustrated in Fig. [T7| We have clarified 
with the authors of Ref. [IB] the numerical origin of the 
discrepancy. This helpful exercise has provided an inde- 
pendent confirmation of the accuracy, within the applied 
physical model, of the theoretical results presented in 
the current pap er. (To compare Fig. [T7] to the original 
plot in Ref. [35| use the conversion to units of normalised 
anisotropy field i?„„/M = 2{E,y^„. - S[iio])/(moM2).) 
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C. Anisotropy fields 

Having analysed the in-plane and out-of-plane 
anisotropies based on the direction of easy axes, we turn 
our attention to the relative strength of the anisotropy 
components, i.e., to the anisotropy energies. The compo- 
nents of magnetocrystalline anisotropy can be described 
in terms of a simple phenomenological model separat- 
ing the free energy density F{M) into components of 
distinct symmetry. Each component is described by a 
periodic function with a corresponding coefficient. We 
find that angular dependencies of the energies obtained 
from our microscopic modelling can be approximated ac- 
curately even in the first and second order of expansion 
into periodic functions of uniaxial and cubic symmetry, 
respectively. 

The coefficients can be determined experimentally, 
e.g., by analysing the FMR spectra i'^^'^^i^^'^^ from 
AMR,^!^ or by fitting SQUID magnetometry data to 
an appropriate phenomenological formula for anisotropy 
energyj ^^i^" In this subsection we extract the relevant co- 
efficients from the calculated anisotropies, track their de- 
pendence on material parameters and compare theory to 
experiment on this level. 

We start with identifying the types of anisotropy terms 
considered in our expansion of the anisotropy energy. 
The cubic anisotropy due to the crystal symmetry of the 
zinc-blende structure is described using terms invariant 
under permutation of the coordinate indices x, y, and 
z. The independent first, second and third order cubic 
terms read: {n^n-y + n^nl + nluy), Kc2 {n^.n^nl), 
and ifc3 {n^^y -f n!^^^ -I- nyufj , respectively, where Ux = 
coscjjsmO, Uy = sin sin 0, and = cos 9 are compo- 
nents of the magnetisation unit vector M (the angles 9 
and (j) are measured from the [001] and [100] axis, respec- 
tively). See the Appendix [Cl for details on the mutual 
independence of all cubic terms. 

As mentioned in previous sections, the cubic 
anisotropy of the host crystal lattice is accompanied by 
different types of uniaxial anisotropy. A generic term 
corresponding to uniaxial anisotropy along a given unit 
vector U depends on the even powers of the dot prod- 
uct (M ■ U). The first and second order terms read: 
Kui{M ■ ilf and Ku2(M ■ C/)"*. The particular cases 
of uniaxial anisotropy terms and their correspondence to 
lattice strains will be described later in this section. 

Before we present the calculated values of the cu- 
bic anisotropy coefficient, we introduce the so called 
anisotropy fields which are often used in literature in- 
stead of the energy coefficients. In this section we plot 
the anisotropy fields in Oersteds (Oe) to make the com- 
parison with experiment more convenient. The relation 
of the anisotropy fields Ha to the energy coefficients Ka 
reads: Ha = 2Ka/M. 

Fig- HH shows Hci and Hc2 as functions of hole density 
p and Mn local moment concentration x at zero tem- 
perature. Both coefficients oscillate as function of the 



hole density p. As discussed in detail in Ref. ^ the 
anisotropies tend to weaken with increasing population of 
higher bands which give competing contributions. Con- 
sistent with this trend the amplitude of the oscillations 
increases with increasing x and decreasing p. The up- 
per limit of the hole density p — Nmu corresponds to no 
charge compensation (Recall, Nmu ~ 2.21a; in 10'^"cm~'^ 
for X in percent). 



2000 n 




FIG. 18: (color online) Lowest order cubic anisotropy field 
Hci and second order cubic anisotropy field Hc2 calculated 
as functions of hole density p (up to zero compensation 
p = Nm„) and Mn local moment concentration x at zero 
temperature; 

Our modelling predicts the extremal magnitude of the 
second order cubic term Hc2 a factor of two smaller than 
the extremal magnitude of the first order term Hd ■ Upon 
increasing the hole density the amplitude of oscillations 
of Hc2 decreases faster than in case of Hd- The third 
order cubic anisotropy field Hcs is negligible compared 
to Hci and Hc2 for all studied combinations of the ma- 
terial parameters. To our knowledge, iJc2 and -ffcs have 
not been resolved experimentally. We emphasise that 
the second order cubic term does not contribute to the 
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anisotropy energy for magnetisation vectors not belong- 
ing to the main crystal plains. The dependence of all 
three calculated cubic terms on the lattice strains of typ- 
ical magnitudes (up to 1%) is negligible. 

Now, we focus on classification of distinct uniax- 
ial anisotropy components and their relation to lattice 
strains lowering the underlying cubic symmetry of the 
zinc-blende structure. We have already mentioned that 
typically the strongest symmetry breaking mechanism is 
the growth strain (introduced in Sec.[n|. It is relevant 
for the in-plane versus out-of-plane alignment of the mag- 
netic easy axis. We have also mentioned the in-plane uni- 
axial anisotropy between the [110] and the [110] axes. Its 
origin is not known, however, we have modelled it using 
the shear strain which is about a factor of ten weaker 
than the typical growth strain. 

Some (Ga,Mn)As epilayers^^iiS, also show a very 
weakly broken symmetry between the main crystal axes 
[100] and [010]. We will introduce here a uniaxial strain 
that can account for this type of anisotropy, however, our 
main motivation for introducing this third strain tensor 
is to complete an in-plane strain basis. This basis is used 
in Sec. llVl to describe all types of lattice in-plane strains 
induced experimentally by growth and post-growth pro- 
cessing of the (Ga,Mn)As epilayers. Once the strain ten- 
sors and corresponding anisotropy contributions to the 
free energy are introduced, it will be shown that the cho- 
sen basis has the advantage of coUinearity of the strain 
and of the resulting anisotropy component. Finally, in 
this subsection the numerical data and comparison with 
experiment will be presented for the bare unpattcrned 
epilayers. The patterned structures will be discussed in 
SecHVl 

Firstly, we recall the growth strain introduced in 
Eq. ([3]). It is usually referred to as the biaxial pseu- 
domorphic strain as it is due to the lattice missmatch 
between the substrate and the epilayer. The doped crys- 
tal is forced to certain dimensions by the substrate in 
the two in-plane directions whereas it can relax in the 
perpendicular-to-plane direction keeping the requirement 
of zero net force acting on the crystal: = Ci2exx + 
Ci2eyy -l- ciiCzz- The corresponding strain tensor: 



resented by a tensor: 
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cii "J 



(9) 



describes an expansion (contraction) along the [100] and 
[010] axes for positive (negative) eo accompanied by a 
contraction (expansion) along the [001] axis. Parameters 
cii and Ci2 are the elastic moduli. The growth strain 
enters our model via the strain Hamiltonian Tistr (see 
Eq. ([T])) and induces a uniaxial anisotropy component 
which can be described in the lowest order by an energy 
term -K[ooi]nl = -if[ooi] cos^ 0. 

The shear strain, first introduced in Sec. Ill Al is rep- 



e = 




(10) 



Positive (negative) k corresponds to turning a square into 
a diamond with the longer (shorter) diagonal along the 
[110] axis. We have used this type of strain as the "in- 
trinsic " shear strain e™* to model the difference in energy 
for magnetisation aligned with the two in-plane diago- 
nals. It results in uniaxial anisotropy along the diago- 
nals, described in analogy to the growth strain by a term 
-~K[iio]{ny-nx f/2 = -K[iio] sin^ (</. - 7r/4) sin^O. 

Finally, we write down the third element of the in-plane 
strain basis: 



e = 




(11) 



Positive (negative) A corresponds to turning a square into 
a rectangle where expansion (contraction) along the [100] 
axis is accompanied by a contraction (expansion) along 
the [010] axis of the same magnitude. Much like in case of 
the growth strain and the shear strain, the requirement 
of zero net force acting on the crystal is kept but this 
time it results in ezz = 0. The strain e" induces uniaxial 
anisotropy along the main crystal axes, described by a 
term — i^[ioo]"y ~ ^-^^^[100] sin^ (/)sin^ 9. 



Let us remark that strain tensors in Eqs. ((9] fTT|l are ex- 
pressed in Cartesian coordinates fixed to the main crys- 
tallographic axes. Strains e** and e" for k = A are related 
by a rotation about the [001] axis by 7r/4, however, the 
cubic crystal is not invariant under such rotation so the 
two strains induce anisotropies with magnitudes i^[ioo] 
and K^iiQ] which are different in general. The growth 
strain e^, the shear strain e**, and the uniaxial strain e" 
can be characterised by a single direction of deformation 
and induce uniaxial anisotropy components aligned with 
that particular direction. We found that higher order 
uniaxial terms are small unless we approach experimen- 
tally unrealistic large values of exchange splitting (large 
x) and hole compensation (low p) . 

In total, we can write our phenomenological formula 
approximating accurately the calculated free energy den- 
sity of an originally cubic system subject to three types 
of strain as a sum of distinct anisotropy components: 



F{M) = {nlnl + nlnl + nlnl) + K,2 (nlnlnl) - 

1x)^-K,mI (12) 



2 -^"^[110] , 

K[mi\nz T^yn, 



By definition of the terms, a positive coefficient -ftr[ooi] 
prefers perpendicular-to-plane easy axis (PEA); positive 
i^liio] and if[ioo] prefer easy axis lying in-plane (lEA) 
aligned closer to [110] and [010] axis, respectively. Note 
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that the anisotropy terms entering the phenomenologi- 
cal formula follow a sign convention consistent with with 
existing literature.— '^'^'^ 

We now provide the microscopic justification for the 
choice of the elements e** and e" of the in-plane strain ba- 
sis and corresponding phenomenological uniaxial terms. 
This will be based on symmetries of the Kohn-Luttinger 
Hamiltonian TLkl and the strain Hamiltonian TLstr as 
shown in Eqs. (|A11|) and (|B2|) . respectively, which re- 
lates the band structure to a general in-plane strain with 
the components Cxx, e-yy, and e^y 

First let us point out that the basis element (the 
growth strain) is invariant under rotation about the [001] 
axis and according to our calculation does not influ- 
ence the in-plane direction of the easy axis (in the linear 
regime of small deformations). We continue by show- 
ing that for and e", the strains and the correspond- 
ing magnetocrystalline anisotropy components are indeed 
coUinear and that this coUinearity applies only for the 
special cases of uniaxial symmetries along the in-plane 
diagonals or main axes. Let us assume a rotation of the 
tensor e" by an arbitrary angle uj about the [001] axis: 

/a o\ 

e"(w) ^ Rl -A i?^ (13) 
y / 

A cos 2u) A sin 2lo 
Asin2cj — Acos2cj 


where is the rotation matrix. (The same analysis ap- 
plies to a rotation of e*). The parameters exx = ~Gyy = 
A cos 2cj and Cxy = A sin 2uj enter the strain Hamiltonian 
(see Eq. (IB2p in the Appendix) only via the matrix ele- 
ment: 



2m 



[72 cos 2(j) — 173 sin 2q 



(15) 



a2 



\/3( 



= —A 02^3 cos(2uj) ~ ia^ sm(2uj) , (14) 

where a2\/3 ^ 03 are strain Luttinger constants. More- 
over, the strain component exy quantifying the shear 
strain enters only Im(c*), whereas the components Cxx ~ 
—Cyy enter only Re(c''). According to our calculation 
the imaginary and real part of c'* generate independent 
uniaxial anisotropy components along the [110] and [100] 
axis, respectively. Their combined effect can be under- 
stood based on an analogy of the in-plane rotation of the 
strain tensor e" and an in-plane rotation of a fc-vector. 

As mentioned in Sec. the Kohn-Luttinger Hamilto- 
nian Ti.K L and the strain Hamiltonian Tistr have the same 
structure. We write here explicitly the matrix component 
c of the Hamiltonian Ti.KL analogous to as a function of 
the in-plane angle of the /c-vector k — |k|[coS(/), sini/), 0]. 
The element reads: 



where again 72 ^ 73 are Luttinger constants describing 
a cubic crystal. For 72 = 73 the Hamiltonian Hkl has 
spherical symmetry. Similarly, if a2\/3 = 03, the strain 
Hamiltonian Tistr is spherically symmetric and the con- 
tributions of Im(c'*) and Re(c'') to the anisotropy of the 
system combine in such a way that the resulting uniax- 
ial term is collinear with the strain e"(u;) rotated with 
respect to the crystallographic axes by an arbitrary in- 
plane angle u. 

Clearly, the underlying cubic symmetry of the host 
crystal causes a non-collinearity of the uniaxial strain 
along a general in-plane direction and the corresponding 
anisotropy component. Moreover, the misalignment is a 
function of Mn local moment concentration, hole density 
and temperature. We discuss further this misalignment 
in more detail in Sec. IIVI Here we point out the distinct 
exception when lu is an integer multiple of 7r/4 and ei- 
ther the real or the imaginary part of vanish rendering 
the strain Hamiltonian effectively spherically symmetric. 
We choose quite naturally the simple forms of e"(ti;) with 
uj = and a; = 7r/4 as elements of the in-plane strain ba- 
sis. For a different choice of the basis elements than in 
Eqs. (fTO|) and (jlip. setting up the phenomenological for- 
mula would be more complicated. 

We can now resume our discussion of the interplay of 
the cubic and uniaxial anisotropy components. Adding 
the uniaxial terms leads to rotation or imbalance of the 
original (cubic) easy axes as shown in Sec. lIIIBl in Fig.[7l 

Fig. [TH shows i7[iio] = 2K[iiq]/M and i7[ioo] = 
2K[iQQ]/M as functions of hole density p and Mn lo- 
cal moment concentration x at zero temperature. Both 
anisotropy fields denpend on material parameters in a 
qualitatively very similar manner. Moreover, we observe 
similar dependence on the doping parameters also in case 
of the field i?[ooi] (not plotted). All the three fields os- 
cillate as functions of hole density. The period of the 
oscillation is longer than in case of Hd- In general, the 
amplitude of the oscillations decreases with decreasing 
Mn local moment concentration. 

The uniaxial fields are linearly dependent on the strain 
from which they originate, unless the strains are very 
large (> 1%). For the shear strain of the value exy — k Ki 
0.01%, which is the typical magnitude in our modelling, 
and zero temperature, the extremal values of -ff[iio] are 
an order of magnitude smaller than the extremal values 
of Hci ~ 10'^ Oe. For typical compressive growth strain 
Eq « —0.2% of an as-grown 5% Mn doped epilayer and 
zero temperature the extremal values of -ff[ooi] are of the 
same order as Hd ■ When the magnitude of the uniaxial 

or 



strain along [100] axis is set to {e^ 



)/2 



--xy, 



V3h^ 

2™ 



b2{kl-kl)-2i{j3kxky)] 



equivalently k = A, -ff[ioo] is approximately a factor of 
two smaller than iJ^xio]- 

To quantify the observed similarity in the calculated 
dependencies of the uniaxial anisotropy coefficients on x, 
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FIG. 19: (color online) Calculated anisotropy fields -ff[iio] and 
i/[ioo] as functions of hole density p (up to zero compensation 
p = Nmv.) and Mn local moment concentration x at zero 
temperature and eo ~ —0.2%. For -^[110] the in-plane strains 
are k = 0.01% and A = {e^y — 0.01%, e^x — Syy — eo), 
while H[iooi ^® found for n — and A — 0.01% (e^y = 0, 
= eo + 0.01%, eyy = eo - 0.01%). 

p, and strains, we can write approximate relationships: 

-^[ooi](a;,P,eo) ~ '?[ooi](a;,p)eo, 
Kiioo]{x,p,X) ~ q[ioo]{x,p)X, 
K[iio]{x,p,k) ~ q[iio]{x,p)K. (16) 

Note, that each anisotropy component depends only on 
one type of strain, which is due to the choice of the basis 
in the strain space (see Eqs. ([9]), (fTO)) . and ifTTjl). (Such 
exclusive dependence of a particular uniaxial anisotropy 
component on the corresponding strain is, indeed, ob- 
tained also from simulations of systems subject to com- 
binations of all three types of strain.) The linearity of 
anisotropy coefficients as functions of lattice strains is 
limited to small elastic deformations of the lattice. The 
approximation cannot be used for strains greater than 1% 



as revealed also by calculations in Ref. [l|. Experiment 
confirms the linear behaviour in case of the growth strain 
up to eo ~ ±0.3%.^- Linear dependence on in-plane uni- 
axial strains is corroborated by experiments discussed in 
Sec.IVl 

In addition to the linearity with respect to strain, 
we observe universal dependence of the three uniax- 
ial anisotropy coefficients on hole density and Mn local 
moment concentration. It can be expressed using the 
anisotropy functions: 

Qlooi]{x,p) ~ q[wo]{x,p) ~ 0A3q[iio]{x,p). (17) 

The anisotropy function q[iiQ]{x,p) due to shear strain is 
approximately twice as large as the anisotropy functions 
q[iOQ]{x,p) and q[ooi]{x,p). A general property of these 
functions is that at medium hole densities a relative com- 
pression yields a tendency of the easy axis to align with 
that direction. On the other hand, for very low and high 
hole densities, the magnetisation prefers alignment par- 
allel to the direction of lattice expansion. 

We caution that Eqs. and pT|) are included to 
promote the general understanding of the anisotropic be- 
haviour of the strained crystal but are not precise. The 
relative error of the approximation given by Eq. (|17p av- 
eraged over the x — p space shown in Fig. is less than 
20%, however, the relative error can be much larger at a 
given combination of x and p where the anisotropy coef- 
ficients fall to zero. 

To finish the analysis of the theoretical results we in- 
clude Fig. [20] to improve the legibility of the data. The 
individual curves correspond to cuts through the 3D 
plots in Figs. [TH] and [H] at fixed Mn local moment con- 
centrations. As already mentioned, the dependence of 
anisotropy fields on hole density is oscillatory. Note that 
the critical hole densities, where the sign inversion oc- 
curs, shift away from the extremal values, i.e., zero hole 
density and zero compensation p = Nmn, with increasing 

X. 

Neglecting the complexity of the dependence of the 
band structure on M (whether changed by doping or 
temperature), one would expect the cubic anisotropy 
coefficient Kd to be proportional to and uniaxial 
anisotropy coefficients -fC[ooi], ^[100] i ^-^d i^inoj to M^. 
In Fig. [201 we can identify intervals of hole density where 
any change in Mn concentration, and therefore in M, does 
not induce a sign change of the anisotropy fields and the 
functional forms of Ka{M) are roughly consistent with 
the above expectations. For other hole density intervals, 
however, the behaviour is highly non-trivial and the func- 
tion Ka{M) can even change sign. 

We now proceed to the discussion of how the theoreti- 
cally expected phenomenology detailed above is reflected 
in experiments in bare unpatterned (Ga,Mn)As epilay- 
ers. The experimental results^'^i^'S are often anal- 
ysed using the following version of the phenomenological 
formula: 

F{M) = -2ttM^ sin^ 9 - K2± cos^ 9 - ^K^^ cos"* 9 
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FIG. 20: Anisotropy fields Hci, -ff[iioij and -ff[ioo] as func- 
tion of hole density (up to zero compensation p = NMn) at 
four Mn local moment concentrations x, zero temperature and 
growth strain eo = —0.2%. For -ff[iio] the in-plane strains are 
K. = 0.01% and A = 0, while -ff[ioo] is found for k = and 
A = 0.01%. (The field Hci is not a function of lattice strains 
but the same values as for calculation of -ff[iio] were used.) 
"Critical" hole densities, where the anisotropy fields change 
sign, are dependent on Mn local moment concentration. 
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where angle 6 and </> are measured, as above, from the 
[001] and [100] axis, respectively. The first term in 
Eq. corresponds to the shape anisotropy described 
in Sec. Ill Bl and not included in Eq. (fT2|) . The uniax- 
ial anisotropy coefficients K2±_ and K2\\ correspond to 
the coefficients JiTjQO]^] and if[iio] in the phenomenologi- 
cal formula Eq. (fT^ . respectively. To identify the third 
and fourth term in Eq. (|18p we rewrite those terms as 
(see also Eq. ICip : 



- / 3 + cos 4(/) . 4 

--Ka\\ sm I 

2 V 4 
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i(if4±-i^4||)cOs4( 

1 



COS 



- (19) 
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where c is an angle independent constant. From here 
we see that the coefficient K^^^^ corresponds to the low- 
est order cubic coefficient Kd in Eq. (|12p and K4± — 
= -fi'[ooi]2 corresponds to the second order uniax- 
ial anisotropy coefficient Ku2 for U \\ [001]. We point 
out that omission of the second order cubic term (and 
other higher order terms) can make the determination of 
^[QOi]o from fitting the data to the phenomenological for- 



mula in Eq. (|T8|) unreliable. Moreover, the accurate ex- 



traction of the coefficient Ki 



[001]: 



can be difficult in sam- 

52 



pies with large value of the first order coefficient X^ooi 
We therefore only note that K^qq^^ extracted from the 
experiment^i^iS never dominates the anisotropy, con- 
sistent with our calculations, and do not discuss the co- 
efficient further in more detail. 

The predicted strong dependence of -^"[001]; ^[110] ? ^^^d 
Kci on hole density, Mn local moment concentration and 
temperature is consistently observed in many experimen- 
tal papers. We start with experiments where the out- 
of-plane anisotropy is studied. Measurements focusing 
mainly on the in-plane anisotropics are discussed at the 
end of this section and in Sec. llVl for patterned or piezo- 
strained samples. 

The coefficient A'[qoi] is extracted in Ref. 58] using 
detailed angle-resolved magnetotransport measurements 
at 4 K for different growth strains in as-grown and an- 
nealed, 180 nm thick samples with identical nominal 
Mn concentration x « 5%. The growth strain rang- 
ing from eo — —0.22% (compressive) to eo — 0.34% 
(tensile) is achieved by MBE growth of (Ga,Mn)As on 
(In,Ga)As/GaAs templates. The observed linear depen- 
dence of -ftTjooi] on Co agrees on the large range of eo with 
the prediction given in Eq. The calculated and 
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measured gradients are of the same order of magnitude 
and sign, and depend on the hole density. The off-set at 
zero strain in the measured dependence of X[ooi] on cq 
in Ref. [Ill is due to the shape anisotropy. 

Ref. [ill presents 50 nm thick, annealed samples with 
nominal Mn doping x = 7%. All the samples are first 
passivated by hydrogen and then depassivated for differ- 
ent times to achieve different hole densities while keeping 
the growth strain the same. The FMR spectroscopy is 
carried out for in-plane and out-of-plane configurations. 
There is qualitative agreement of calculation and mea- 
surement on the level of the directions of the easy axes as 
discussed in the previous subsection. The sign change of 
the uniaxial anisotropy fields driven by increase of tem- 
perature is observed. The measured coefficients -^"[001] 
and Kci are of the same order of magnitude as the calcu- 
lated ones and ^^[ooi] ~ K^i is consistent with the weaker 
growth strain in annealed samples. 

Ref. [stI ] presents an as-grown, 6 nm thick film nomi- 
nally doped with Mn to x = 6%, grown on Gao.76Alo.24As 
barrier doped with Be. Increasing the Be doping in- 
creases the hole density without changing the Mn local 
moment concentration. The fitting of the FMR spectra 
is done using the coefficients -ft'iooi] a-nd Kd and the g- 
factor of the Mn. The anisotropy field corresponding to 
the coefficient K^qq^ reaches value as high as « 6000 Oe 
at 4 K. Large values of -fir[ooi] is consistent with expected 
large growth strain in a thin as-grown samplei^i^ How- 
ever, for the measured -ftr[ooi] our calculations would im- 
ply strain eo ~ 1% which is an order of magnitude larger 
than typical strains in as-grown x = 6% (Ga,Mn)As ma- 
terials. Other effects are therefore likely to contribute 
to ^^[001] in this sample. (Confinement effect or inhomo- 
geneities are among the likely candidates.) The experi- 
mental if [001] (^ci) increases (decreases) with increasing 
hole density which is in agreement with our modelling of 
highly compensated samples. 

Observation of qualitatively consistent behaviour of 
the anisotropics with the theory but unexpectedly large 
magnitudes of the anisotropy fields appli es also to thick 
samples studied by FMR in Refs. [53156]. Temperature 
dependence of the anisotropy fields is studied by FMR 
in Ref. [56| for a low doped {x « 2%), as-grown, 200 nm 
thick (Ga,Mn)As film. Only the combined contribution 
of shape anisotropy and K[qqi] was resolved. The easy 
axis stays in-plane for all studied temperatures which is 
consistent with predicted crystalline anisotropy as well as 
the shape anisotropy dominating at weak growth strains. 
The uniaxial in-plane anisotropy is of the predicted mag- 
nitude but its sign corresponds to modelling by the less 
frequent negative intrinsic shear strain. 

Ref. [131 discussed in Sec. IIII Al on the level of easy 
axis orientation shows, among other samples, 300 nm 
thick annealed epilayers with nominal Mn concentration 
x = 3% deposited on GaAs and (Ga,In)As substrate un- 
der compressive and tensile growth strain, respectively. 
The strain is measured by x-ray diffraction, however, 
the predicted linear dependence of i^iooi] on the growth 



strain (Eq. ((T6)) ) cannot be tested due to different sat- 
uration mag netisation and Tc in both samples. Both 
Refs. [S^IsgI I report the coefficient Kd in the 300 nm and 
200 nm thick samples an order of magnitude larger than 
the calculated one which can-^ be attributed to sample 
inhomogeneities in these thick epilayers. Ref. [H] stud- 
ies also 120 nm thick, annealed and as-grown epilayers 
with X = 8% deposited on GaAs. The coefficient if[ooi] 
doubles its value at low temperature on annealing. Both 
-/^[ooi] and Kd in the thinner samples have values of the 
order predicted by theory for material with Mn doping 
X = 8%. 




FIG. 21: (Color online) Angle tp of the easy axis with respect 
to the [110] axis as function of hole density p (up to zero 
compensation p = Nmti) and Mn local moment concentration 
X at zero temperature, eo = —0.2%, and k — 0.01%; 

Now we analyse experiments focusing on the in-plane 
anisotropy where the relevant anisotropy coefficients are 
Kd and ifjuo]. Note that the experimental papers dis- 
cussed below mostlyiii^i^i^i2£ use the notation with the 
in-plane magnetisation angle tp measured from the [110] 
axis. To avoid any confusion we write the in-plane form 
of Eq. p2| using the original anisotropy coefficients and 
the angle tp = (p + tt/A: 

F{M) = - ^ sin^ 2iP + K^^q^ sit? ^. (20) 

To facilitate the comparison with experiment we use the 
notation of Eq. ()20|) consistently in the remaining parts 
of this paper. 

The magnetic easy axes lie closer to the [100] or 
[010] direction than to any diagonal when Kd > and 
V2-^[iio] < ^ci- Negative Kd always leads to diagonal 
easy axes. We include Fig. [21] to elucidate the combined 
effect of Kd and i^ino] on the in-plane direction of the 
easy axes. The angle tpEAix,p), plotted as a function 
of Mn concentration and hole density at zero tempera- 
ture minimises the free energy F{M). The local minima 
at -0 = 0° (black) and ip = 90° (white) are formed for 
negative Kd- When K^n^^ is positive (negative), the 
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global minimum is at ^/j = 0° {ip = 90°). The higher 
energy local minimum disappears for |i^ci| = l^[iio]l- 
Only one energy minimum forms for \Kci\ < |i^[iio]| and 
for positive (negative) if [no] the easy axis is at V' = 0° 
{ip = 90°). The interface of black and white regions is an 
evidence of a discontinuity of the function ipEAiXjp) due 
to switching of the sign of K^hq] when Kd < 0. The grey 
(coloured online) regions in Fig. [^T] correspond to com- 
petition of cubic and uniaxial anisotropy when Kd > 
and \Kci\ > |if[iio]|. Then there are two easy axes at 
^Pea and 180° — ipEA forming "scissors" closing at the 
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FIG. 22: Calculated anisotropy fields K^i and if[iio] as func- 
tion of temperature and magnetisation at two hole densities 
(given in units u = 10^° cm"''), Mn concentration x — 2%, 
strains eo = —0.2%, n — 0.005%. Irregular behaviour is ob- 
served for the lower hole density. 

To demonstrate the typical scaling of in-plane 
anisotropy components with temperature, we discuss the 
50 nm thick as-grown (Ga,Mn)As epilayer with Mn con- 
centration X — 2.2% determined by x-ray diffraction and 
secondary ion mass spectrometry, presented in Ref. [36| . 
The anisotropy coefficients if [no] and Kd are obtained 
by fitting to the M{H) loop with magnetic field along 
the hard direction. They can be compared to Fig. [221 



which shows the calculated anisotropy fields as func- 
tions of temperature for two values from the interval of 
hole densities corresponding to the as-grown sample. For 
p — 2.5 X 10^° cm~'^ both the calculated and measured 
Kd is greater than if [no] at low temperatures but be- 
comes smaller than if [no] at T w Tc- The calculated 
ifci is an order of magnitude smaller than the experi- 
mental one, however, there is agreement on the level of 
the temperature dependent ratio of ifd and if [no]. On 
the contrary. Fig. [22] shows a non-monotonous depen- 
dence of ifci on temperature for p = 1.5 x 10^° cm~'^. 
This singular behaviour is not measured in Ref. 13611 but 
it is reported in a more systematic study in Ref. [4l|. 

The temperature dependence of anisotropy coefficients 
if [no] and ifd is studied by planar Hall effect in Ref. [6^. 
The mutual behaviour of the two coefficients observed 
in the as-grown (Ga,Mn)As epilayer with nominal Mn 
concentration x « 4% and Tc = 62 K is qualitatively the 
same as in Ref. [36j . ifd becomes smaller than if [no] 
at T = 26 K which is in agreement with our modelling. 
No sign change of ifd is reported in this experimental 
work. Again, the calculated ifd is an order of magnitude 
smaller than the experimental one. 

Ref. jlJl resolves the in-plane coefficients ifd and 
if[iio] in four samples with nominal Mn doping x = 7% 
and different hole densities. In samples with lower hole 
densities the dependence of ifd and if [no] is qualita- 
tively consistent with Ref. [s^, however, both coefficients 
change sign when temperature is increased in samples 
with higher hole densities (p 10^^ cm~^, Tc = 130 K). 
Our model predicts such sign change for a short interval 
of high hole compensations and a larger interval of low 
hole compensations as shown in Fig. I^Of a). 

Another type of temperature scaling of ifd and if [no] 
is observed in a 50 nm thick, annealed sample with nom- 
inal Mn doping X = 7% and Tc = 165 K.^° iC[no] 
is larger than ifd on the whole temperature interval 
{T — 4 ~ 165 K). Both coefficients are positive, de- 
crease on increasing temperature, and their magnitudes 
are of the same order of magnitude as the calculated 
anisotropics. The stability of sign of if [no] is observed 
theoretically for higher "intrinsic" shear strain as dis- 
cussed in Fig. [Tljin Sec. IIIIBI 

The temperature dependence of domain wall proper- 
ties of a 500 nm, as-grown (Ga,Mn)As ffim with Mn dop- 
ing X — 4% is studied by means of the electron hologra- 
phy in Ref. |17| • The width and angle of the domain walls 
were determined directly from the high-resolution im- 
ages. The ratio of the anisotropy coefficients if [no] /ifd 
was extracted from these observations combined with 
Landau-Lifshitz-Gilbert simulations. The Neel type do- 
main walls evolve from near-90°-walls at low tempera- 
tures {T — 10 K) to large angle [110]-oriented walls and 
small angle [110]-oriented walls at higher temperatures 
(T = 30 K). The angles of domain walls aligned with 
particular crystallographic directions reveal positions of 
the magnetic easy axes. The "scissors" of the easy axes 
(described in discussion of Fig. [21]) are closing around the 



20 



[110] axis on increasing temperature consistent with our 
modelling. 

The domain- wall width is inversely proportional to the 
effective anisotropy energy barrier between the bistable 
states on respective sides of the domain wall: K'^(iff^ = 
Kci/'i - K{iw]/2 ([110]-oriented walls) and K^^^J^^ = 
+ iv:[iio]/2 ([ITO] -oriented walls). The width of 
the [110]-oriented wall in Ref. [l3| initially increases 
with temperature and then saturates at high tempera- 
ture while the [110]-oriented wall width keeps increasing 
with temperature until it becomes unresolvable. This 
observation corresponds well to the theoretical predic- 
tion and can be qualitatively understood by consider- 
ing the approximate magnetisation scaling of Kci ^ M^, 
if [110] ~ M^, and magnetic stiffness ~ Af^. 

Finally, Refs. [Hm] present (Ga,Mn)As field-effect 
transistors (FETs), where hole depletion/accumulation 
is achieved by gating induced changes of the in-plane 
easy axis alignment. In Ref. [5§| the Mn doped layer 
is 5 nm thick with Mn doping x = 2.5% and hole density 
p ~ 1 X 10^^ — lO^^cm"^. The direction of magnetic easy 
axes was detected by AMR at T = 4K. The 20% variation 
of the hole density achieved by applying the gate voltage 
from —IV to 3V is determined from variation of the chan- 
nel resistance near Tc- This value was a starting point 
for simulations of the depletion at T = 4 K giving hole 
density changes Ap w 5 x lO^^cm^'^. The measured Kd 
is negative and its magnitude decreases with depletion. 
The theoretical magnitude (~ 10 mT) and sign of Kd 
for the relevant hole density range, as well as the varia- 
tion of Kci with varying hole density, are consistent with 
the experiment. Recall that negative Kd corresponds to 
diagonal easy axes captured by two black /white regions 
in Fig. H H Sa mples reported in Refs. [ioillliiil (see 
also Sec. IIII B[) and in Ref. [10] with diagonal easy axes 
at low temperatures fall into the right region with lower 
hole compensations, whereas the sample in Ref. [H^ is 
a rarely observed example of diagonal easy axes at high 
compensation and low temperature corresponding to the 
left black/white region in Fig. [21] 



IV. SAIVIPLES WITH POST-GROWTH 
CONTROLLED STRAINS 

In the previous section, we discussed three types of 
lattice strain and calculated corresponding types of uni- 
axial anisotropy components. In the bare, unpatterned 
epilayers we could analyse and compare to experiment 
only anisotropics induced by the growth strain and by 
the unknown symmetry breaking mechanism modelled 
by the "intrinsic" shear strain. The calculations includ- 
ing the model shear strain allow us also to estimate the 
magnitude of real in-plane lattice strains, controlled post- 
growth by patterning or piezo stressing, that can induce 
sizable changes of anisotropy. In this section we inves- 
tigate samples where these post-growth techniques are 



used to apply additional stress along any in-plane direc- 
tion. We will focus primarily on stresses along the main 
crystal axes and in-plane diagonals. We will also com- 
ment on the procedure for determining the lattice strain 
from specific geometrical parameters of the experimental 
setup. Where necessary, we distinguish the externally in- 
duced strain and the "intrinsic" shear strain, which mod- 
els the in-plane symmetry breaking mechanism already 
present in the bare epilayers. Returning to the nota- 
tion of Sec. |TT] we denote the latter strain by the symbol 
e™* . For better physical insight and to relate with discus- 
sion in previous section we will map the anisotropics on 
the phenomenological formulae by decomposing the total 
strain matrix into the three basis strains (Eqs. ([ MTT]) '). 
We will then write the corresponding anisotropy energy 
terms as in Sec. IIII CI assuming linearity between the re- 
spective basis strains and anisotropy energy components 
(see Eq. (UHl)). Experiments will be discussed based on 
microscopic anisotropy calculations with the total strain 
tensor directly included into the Hamiltonian. 

We begin this section by discussion of the in-plane uni- 
axial strain induced by post-growth lithography treat- 
ment of Mn-doped epilayers grown under compressive 
lattice strain. Narrow bars with their width comparable 
to the epilayer thickness allow for anisotropic relaxation 
of the lattice matching strain present in the unpatterned 
film. An expansion of the crystal lattice along the di- 
rection perpendicular to the bar occurs while the epi- 
layer lattice constant along the bar remains unchanged. 
Parameters sufficient for determination of the induced 
strain are the initial growth strain eo and the thickness 
to width ratio t/w of the bar. In the regime of small 
deformations the components of the induced strain are 
linearly proportional to the growth strain. The strain 
tensor for a bar oriented along the [100] axis reads: 

/ + 1 \ 
efioo] = eo 1 , (21) 

V ofi7(^-2); 

where the lattice relaxation is quantified by p which is a 
function of t/w and can vary over the bar cross-section. 
We calculate the distribution of p over the cross-section 
of the bar using Structural Mechanics Module of Com- 
sol (standard finite element partial differential equation 
solver, www.comsol.com). Since the macroscopic simu- 
lations ignore the microscopic crystal structure, they ap- 
ply to bars oriented along any crystallographic direction. 
We therefore introduce a coordinate system fixed to the 
bar: a;'-axis lies along the relaxation direction transverse 
to the bar, y'-axis along the bar, and z'-axis along the 
growth direction. We approximate the bar by an infinite 
rectangular prism with translational symmetry along the 
j/'-axis, attached to a thick substrate. 

Fig. [23] shows the spatial dependence of the function 
p{x',z') for a given thickness to width ratio and com- 
pressive growth strain eo < 0. Only the area of the bar 
is plotted, whereas the strain induced in the patterned 
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FIG. 23: Spatial dependence of the strain coefficient p due 
to lattice relaxation in a narrow bar with t/w = 0.4 and 
compressive growth strain eo < 0, simulated values of p are 
plotted for the cross-section of the bar. 
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part of the substrate is not shown. (The substrate re- 
laxation is not directly related to the microscopic simu- 
lation of the anisotropy energy). In wide bars {t/w ^ 1) 
the relaxation is very non-uniform, whereas narrow bars 
{t/w ^ 1) are fully relaxed. Fig. [24] shows still a fairly 
non-uniform relaxation for t/w = 0.4 with large relax- 
ation at the edges. We point out in this case that the 
resulting anisotropy can be very sensitive to the details 
of the etching (vertical under-cut/over-cut profile). 

The non-uniform strain distribution in wider bars can 
in principle force the system to break into magneti- 
cally distinct regions. However, experiments show rather 
that the whole bars behave as one effective magnetic 
medium. Because of the linearity between the strain and 
the anisotropy (see Eq. (dS])) we can model the mean 
magnetic anisotropy by considering the spatial average 
of e^j^Q(,j over the bar cross-section. The inset of Fig. [24l 
shows the averaged value p as a function of the width to 
thickness ratio. It confirms that the effect of relaxation 
can reach magnitudes necessary to generate significant 
changes in the magnetic anisotropy. In very narrow bars 
the induced uniaxial anisotropy can override the intrinsic 
anisotropics of the unpatterned epilayer and determine 
the direction of the easy axis. 

If the bar is aligned with the [100] or [010] crystal 
axis, the strain e|.^ppj in Eq. (PT|) with the average re- 
laxation magnitude p can be used directly as input pa- 
rameter of the microscopic calculation (see Eq. (jB2[) in 
the Appendix) . Alternatively, the total strain tensor can 
be decomposed into the growth basis strain from Eq. ([9| 
and the uniaxial basis strain introduced in Eq. (|11|) : 



3[ioo](eo,p) 
eo 
A 



eS(go) + e"(A), 

P 

-eo-. 



(22) 
(23) 

(24) 



Their effects on the magnetic anisotropy can be consid- 
ered separately utilising the results shown in Sec. IIII CI 
Now we discuss the introduction of uniaxial in-plane 
anisotropics by a piezo actuator attached to the sam- 
ple. In this case, the (Ga,Mn)As film is assumed to 



FIG. 24: (Color online) Sections of p{x' , z) in Fig.[23]at fixed 
values of z (given next to the curves in relative units) of a 
thin bar. Inset shows the average strain 'p{t/w) as a function 
of the thickness to width ratio. 



follow the deformation of the stressor. (The substrate 
is usually thinned to achieve better transmission of the 
piezo-strain to the studied epilayer. Macroscopic Comsol 
simulations predict transmission of approximately 70% of 
the piezo-strain in a substrate with thickness to lateral 
size ratio t/l k, 0.1 and transmission of approximately 
90% of the piezo-strain for t/l « 0.02.) The net effect 
of the piezo-stressing on normal GaAs epilayers has been 
investigated experimentally for example in Ref. [s^ ] for a 
standard PbZrTiOa (PZT) piezo actuator. The induced 
strain can reach magnitudes ~ 10~^ at low temperatures, 
which are sufficient to induce observable anisotropics in 
(Ga,Mn)As, as shown in Sec. IIII CI The deformation is 
linearly proportional to applied voltage on the transducer 
and increases with increasing temperature. 

The dependence of uniaxial anisotropics due to addi- 
tional piezo-strains is analogous to the behaviour of re- 
laxed microbars, however, the form of the strain tensor 
induced by the stressor is typically more complex. Let 
us first assume a strain tensor with components in the 
Cartesian coordinate system fixed to the orientation of 
the piezo stressor: x'-axis lies along the principal elon- 
gation direction, ^'-axis is perpendicular to plane of the 
thin film. We denote the deformation along the x'-axis by 
a and the simultaneous deformation along the y'-axis by 
a' . Note that shear strains are typically not considered 
when describing the action of a piezo-stressor. The third 
parameter describing the strained (Ga,Mn)As epilayer is 
the growth strain eo- Our analysis takes into account 
only structures that can be parametrised by these three 
values. The strain tensor in the dashed coordinate sys- 
tem reads: 



/cr + eo 

cr' + eo 1(25) 

\ -f^(2eo + CT + a' 

Components of this tensor are considered uniform in the 
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studied epilayer. If the principal elongation direction 
of the piezo stressor is aligned with the [100] crystallo- 
graphic axis the strain tensor s^iqq^ can be used directly 
as an input of the microscopic simulation. Similarly to 
the strain induced by lattice relaxation, g^iqq^ can be de- 
composed into the growth basis strain and the uniaxial 
basis strain: 



ePoo](eo,a,a') = eS(go)+e"(A) 
eo 
A 



(26) 
(27) 
(28) 



Again, the results shown in Sec. IIII CI can then be 
used when analysing the resulting magnetocrystalline 
anisotropics. Recall that has a minor effect on the 
in-plane anisotropy and can therefore be omitted when 
discussing in-plane magnetisation transitions. 

So far we have described induced strains aligned with 
the [100] crystal axis. In case of a lattice relaxation or 
piezo stressor aligned at an arbitrary angle tu, the fol- 
lowing transformation of the total strain tensor er,„„, or 



'[100] 



[100] 

to the crystallographic coordinate system applies: 



r(p) 



f/T r(p) p 
^^^'^[100]-^'- 



where the rotation matrix reads: 

cos(cj — 7r/4) sin(a; — 7r/4) 

Rui = I - sin(ti; - 7r/4) cos(lj - 7r/4) 

1 



(29) 



(30) 



The angular shift by — 7r/4 is because we measure the 
angle uj from the [110] axis. This convention was intro- 
duced in Sec. IIII CI before Eq. (f20|) and is used consis- 
tently in this section for all in-plane angles. The rotated 
total induced strain can be used directly as the input 
strain matrix for the microscopic calculation or it can be 
decomposed into all three elements of the in-plane strain 
basis. In case of the relaxation-induced strain, we obtain: 



eo 
A 



-eo^ sin 2lj, 
P 



k = en— cos2ti;. 
2 



(31) 
(32) 

(33) 
(34) 



In case of the rotated piezo stressor, the same decompo- 
sition follows, however, the effective strain magnitudes A 
and k depend on different real experimental parameters: 

eP(eo,fT,a') - eS(go)+e"(A)+e^(«:), (35) 

eo = eo + (36) 



sin 2lij, 



■ cos 2lo. 



(37) 
(38) 



Considering the linear dependence of the anisotropy 
coefficients on the corresponding strain elements (see 
Eq. (|16p ). we can write the part due to post-growth in- 
duced strains of the phenomenological formula for the 
free energy as a function of angles V' and lo: 



(39) 



Fu{M) = if[no](^^)sin^V^-t- 

+/V[ioo](w)sin2(?/' + 7r/4) 

- '?[iio]K(t^)sin^-0 + 

+9[ioo]A(w)sin^(V' + 7r/4), 



where we use the notation analogous to Eq. (|16p in 
Sec. IIII CI The relation of the effective parameters A 
and k to the experimental parameters of microbars or 
stressors oriented along arbitrary crystallographic direc- 
tion is given by Eqs. p3ll34p or p7ll38p . respectively. The 
linearity of the anisotropy constants ifjioo], ^[iio]: £^nd 
/^[QO]^] on corresponding strain coefficients and the form 
of the strain tensors in Eqs. PT|) and ([55]) allow us to fac- 
tor out the cj-dependence of KuS. Figs. [THl [121 and 
together with Eqs. ((3T|) and ([SS]) can therefore be used 
for analysing magnetic anisotropics induced by micropat- 
terning or piezo stressors oriented along any crystallo- 
graphic direction. 

The full angular dependencies of the anisotropy en- 
ergy calculated directly from the total strain tensor in- 
cluded into the Kohn-Luttinger kinetic-exchange Hamil- 
tonian for several combinations of k and A are plotted in 
Fig. [25l Recall that analogous in-plane angular depen- 
dencies of the anisotropy energy were presented in Fig. [5] 
- 1101 where only the competition of the growth strain 
and shear strain e" with the cubic anisotropy of the host 
lattice was considered. 

Fig. [^5f a) shows four angular dependencies of the 
anisotropy energy for x — 3% and p ~ i y. 10^*'cm^'^. 
The curves are marked by the values of the effective 
strain components. The solid curve for weak shear strain 
k = 0.01% and no uniaxial strain A = has two local 
minima close to the main crystal axes indicating domi- 
nant cubic anisotropy with Kd > for the considered 
X and p. The easy axes are shifted due to the positive 
shear strain towards the [110] axis which is the direc- 
tion of relative lattice compression, consistently with the 
discussion in Sec. IIII CI for samples with medium hole 
densities. Additional uniaxial strain A = —0.05% results 
in only one global minimum easy axis rotating towards 
the [100] direction which is again the direction of relative 
lattice compression. 

The dashed curve in Fig. ISST a) corresponding to strong 
shear strain k = 0.09% and no uniaxial strain A = has 
only one global minimum at the [110] diagonal, indicating 
domination of the uniaxial anisotropy over the underlying 
cubic anisotropy. Addition of the uniaxial strain A = 
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(a) a; = 3%, p = 3 X 10^°cm-^ eo = -0.3% when 
A/2 = 0, eo = -0.25% when A/2 = 0.05% 




(b) A = -0.05%, k = 0.01%, eo = -0.25% 

FIG. 25: (Color online) Magnetic anisotropy energy AE — 
-E^— i5[ioo] as a function of the in-plane magnetisation orienta- 
tion M — |M| [cos sin-!/), 0] and its dependence on material 
parameters. Effects of the shear strain and the uniaxial strain 
combine linearly (a). Magnetic easy axes (marked by arrows) 
change their direction upon change of Mn local moment con- 
centration X, and hole density p (in units u = 10'^'' cm"'') for 
a fixed uniaxial and shear strain (b). Both plots assume zero 
temperature. 



—0.05% leads to rotation of the easy axis towards the 
direction of relative compression ([100] for A < 0). 

Curves plotted in Fig. [HJb) differ in the material pa- 
rameters but share the same weak shear strain k — 0.01% 
and the same uniaxial strain A = —0.05%. The solid 
curve for x — 3% and p = A x lO^^cm"^ falls into the 
range of hole densities where the cubic anisotropy coef- 
ficient Kci is positive so the easy axes in the absence of 
in-plane strains align parallel to the main crystal axes. 
Adding the uniaxial strain A yields only one global mini- 
mum along the [100] direction and the shear strain shifts 
the easy axis towards the [110] diagonal. Again, for both 
strains the easy axes tend to align along the direction of 



lattice compression for these medium doping parameters. 

The dashed curve in Fig. I^b) for x = 5% and 
p = 8 X 10^°cm~^ can be described by a negative Kd cor- 
responding to diagonal easy axes in the unstrained bulk 
epilayer. The additional shear strain k makes the [110] 
direction the global minimum easy axis. Note that for 
these values of x and p the easy axis prefers to align with 
the direction of lattice expansion. Consistently, the uni- 
axial strain A rotates the easy axis towards the direction 
of relative lattice expansion, i.e., towards the [010] axis. 
Finally, the dash-dotted curve for x = 5% and high hole 
density p = 12 x 10^°cm~^ corresponds to positive Kd 
and again, when the in-plane strains are included the easy 
axes prefer the direction of relative lattice expansion. To 
summarise the discussion of Figs. I^Sl a) and (b), the pre- 
ferred alignment of the in-plane easy axis with either the 
lattice contraction or expansion direction depends on x 
and p. For a given doping it has always the same sense 
for both the shear strain k and the uniaxial strain A and 
is uncorrelated with the sign of the cubic anisotropy com- 
ponent. These conclusions are independent of the growth 
strain, at least for its typical values eo < 1%. 

Now we analyse experimental studies that control the 
in-plane strain by means of post-growth lithography. 
Refs. [1^ and [l^] present structures with the shear and 
uniaxial strain induced locally by anisotropic relaxation 
of the compressive growth strain. Ref. [13] studies an L- 
shaped channel with arms aligned along the [110] and 
[110] directions patterned by lithography in a 25 nm 
thick (Ga,Mn)As epilayer with nominal Mn concentra- 
tion X — 5%. Hole density p ~ 5 x 10^°cm~^ was es- 
timated from high-field Hall measurements. This pat- 
terning allows relaxation of the growth lattice matching 
strain in direction perpendicular to the channel. There- 
fore, the generated uniaxial strains in each arm of the L- 
shaped channel have opposite signs. The induced shear 
strain is added to (subtracted from) the "intrinsic" shear 
strain in the arm fabricated along the [110] ([110]) axis. 
The magnitude of the induced strain increases with de- 
creasing width of the channel. A large effect on magnetic 
easy axes orientation has been observed in a 1 ^m wide 
channel while only moderate changes have been found in 
a 4 /im bar. In both cases the easy axes of the unpat- 
terned epilayer rotated in the direction perpendicular to 
lattice expansion. The sense and magnitude of the easy- 
axis reorientations in the relaxed microbars are consistent 
with theory prediction for the relevant values of x, p, and 
microbar geometry. 

Refs. [4)] and |26| show lithographically induced uniax- 
ial anisotropy along the [100] or [010] axis in arrays of nar- 
row bars. Ref. [3] presents 200 nm wide bars fabricated 
in an as-grown 70 nm thick film with Mn concentration 
X — 2.5% determined by x-ray diffraction. Ref. [1^ re- 
ports lattice relaxation in 200 nm wide, 20 nm thick bars 
in an as-grown material with nominal Mn concentration 
X = 4%. In both studies the unpatterned epilayers have 
two equivalent easy axes close to main crystal axes. Af- 
ter the anisotropic relaxation of the growth strain in the 



nanobars the easy axis corresponding to the relaxation di- 
rection is lost, whereas the other easy axis is maintained. 
This behaviour is in agreement with our simulations on 
the relevant interval of dopings and patterning induced 
strains. 

The anisotropies induced in the relaxed structures in 
Refs. [i'26''27'' 28i| can be predicted using the results of 
Sec. IIII CI directly. Bearing in mind the negligible effect 
of the growth strain, the relevant part of the strain ten- 
sor describing the relaxation along the main crystal axes 
has the form of the uniaxial basis strain e", as shown in 
Eq. (f22|) . and corresponds to the anisotropy component 
with the previously calculated coefficient K^iqq^ . The re- 
laxation along the diagonals is described by the strain 
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again neglected the contribution from the growth strain FIG. 26: (Color online) Calculated magnetic anisotropy en- 
e^. It induces uniaxial anisotropy component quantified ergy AE = E^- E^^j^^ as a function of the in-plane magneti- 
by the coefficient K\iin\. Note that the "intrinsic" shear sation angle ip measured from the [110] axis at T= 5/8Tc 



strain e™ in the modelling is independent of the exter- 
nally introduced lattice distortion and needs to be added 
to the total strain tensor if the corresponding anisotropy 
is present in the unpatterned epilayer. As mentioned be- 
fore, the simulated rotation of easy axis directions in the 
relaxed microbars is in good agreement with the mea- 
sured behaviour. 

The piezo-strain is also applied in most cases along the 
main crystal axes or diagonals. In Ref. a PZT piezo- 
electric actuator is attached to a 30 nm thick (Ga,Mn)As 
epilayer grown on a GaAs substrate thinned to 100 /ini. 
The principal elongation direction of the actuator is 
aligned with the [110] crystallographic direction. The 
nominal Mn concentration of the as-grown epilayer is 
4.5%. The relative actuator length change is approxi- 
mately 4 X 10^^ at T = 50 K (measured by a strain 
gauge) for the full voltage sweep (from -200 V to 200 V). 
Such piezo-strain induces a rotation of the easy axis by 
/S.ipEA ~ 65°. Our modelling predicts /S.iIjea of the same 
order for relevant material and strain parameters. The 
easy axis rotates towards the [110] ([110]) direction upon 
contraction (elongation) along the [110] axis in agreement 
with the behaviour observed in the relaxed microbars and 
with our modelling. 

Ref. I30l extends the piezo-stressed (Ga,Mn)As study 
in Ref. [2^ to low temperatures. Again, PZT piezo ac- 
tuator is attached to a Hall bar along the [110] crystallo- 
graphic direction. The 30 nm thick, as-grown (Ga,Mn)As 
epilayer grown on GaAs substrate has nominal Mn con- 
centration 4.5% and Tc = 85 K. A strain gauge measure- 
ment shows almost linear dependence of the piezo-strain 
in the Hall bar on temperature (in the range 5 K to 50 K). 
The anisotropy coefficients ^^[iio] and Kd are extracted 
from the angle-dependent magnetoresistance measure- 
ment as a function of temperature for three voltages (- 
200 V, V, and 200 V). At high temperatures the rel- 
ative elongation of the structure is again approximately 
4 X 10~^ and the corresponding uniaxial anisotropy dom- 
inates over the intrinsic uniaxial anisotropy along the 
[110] axis. Close to 5 K the action of the piezo is negli- 



eo = 0, e!,"* = 0.017%, x = 5%, and p = 5 x 10^°cm- 



The 



curves are labelled by o, the induced strain along the princi- 
pal elongation direction of the piezo tilted by angle u — —10°, 
and by the corresponding vol tag e. (The relationship of a and 
voltage is inferred from Ref. [sj] to allow for direct compari- 
son with experiment.) The easy axis rotates smoothly upon 
sweeping the voltage. For -100 V a shallow local energy min- 
imum forms due to the underlying cubic anisotropy (marked 
by arrow). 



gible SO the intrinsic uniaxial anisotropy is stronger than 
the induced one, however, the total in-plane anisotropy 
is dominated by the cubic anisotropy. The measured and 
calculated induced anisotropy along the [110] direction 
are of the same sign and order of magnitude for the con- 
sidered temperatures. 

Ref. [3^ presents a 15 nm thick, annealed sample 
doped to a; = 8%, subject to piezo stressing along the 
[010] axis. The anisotropy coefficients are extracted 
from transverse AMR. The PZT actuator induces rela- 
tive elongation ranging from 1.1 x 10^'^ for voltage 200 V 
to 0.7 X 10^'^ for -200 V, measured by a strain gauge. The 
difference of the limits is again approximately 4 x 10~^ 
but all values are shifted towards tensile strain most likely 
due to different thermal dilatation in the sample and the 
actuator. The lattice expansion along the [010] direc- 
tion leads to alignment of the easy axis along the [100], 
in agreement with our modelling and with the experi- 
mental studies discussed in this section. The extracted 
cubic anisotropy field is roughly a factor of two lower 
compared to studies of samples with high hole compen- 
sation sharing the value « 1000 Oe at different nominal 
Mn concentrationsi^SiSi^ The low critical temperature 
Tc = 80 K suggests lower effective Mn concentration 
in Ref. Our calculations for lower Mn local mo- 

ment concentration and high hole compensation predict 
the anisotropy coefficients Kd and ^^[loo] induced by the 
piezo strain in correspondence with the measured coeffi- 
cients. 

Finally, we discuss a piezo-strain induced along a gen- 
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eral in-plane direction. In Ref. 'Sl'l the principal elon- 
gation direction of the PZT piezo actuator is tilted by 
angle uj = —10° (with respect to the [110] axis). The 
25 nm thick, as-grown (Ga,Mn)As epilayer with nomi- 
nal Mn concentration x — 6% is grown on a GaAs sub- 
strate, which was thinned before attaching of the stres- 
sor to « 150/im. The anisotropics are determined from 
SQUID and AMR measurements at 50 K. The uniaxial 
strain caused by differential thermal contraction of the 
sample and the piezo on cooling (at zero applied voltage) 
is of the order 10"''. The uniaxial strains generated at 
the voltage ±150 V are cr fa ±2 x 10"'' and a' « —a/2 at 
50 K. The tilt of the piezo with respect to the crystal di- 
agonal results in a complicated interplay of the intrinsic 
and induced anisotropy. The easy axis of the bare sample 
aligns with the [110] axis due to strong intrinsic uniaxial 
anisotropy with -ft'ino] > Kd > 0. This easy axis rotates 
to an angle tpEA — 65° upon attaching of the piezo and 
cooling to 50 K. Application of -1-150 V to the stressor 
causes the easy axis to rotate further to ^pE A = 80° while 
for -150 V the axis rotates in the opposite direction to 
'0-EA = 30°. Note that the negative voltage weakens the 
total piezo-strain and allows domination of the intrinsic 
anisotropy with easy axis closer to the [110] axis. 

The hole compensations expected in Ref. [3l| are in 
the range p/Nun = 0.6 — 0.4 and the relevant range of ef- 
fective Mn concentrations is a; = 3 — 5%. i^Tjuo] measured 
in the bare epilayer is modelled by e™* = 3 — 2 x 10^'' 
(slightly weaker than the strain induced in the structure 
at zero piezo- voltage). The in-plane anisotropy ener- 
gies calculated on this parameter interval using the total 
strain tensor (induced and "intrinsic" components) are in 
good quantitative agreement with the easy axis orienta- 
tions measured at the three piezo voltages. Fig. shows 
calculated curves for one representative combination of x, 
p, and e™* from the relevant interval, for the fixed tilt of 
the stressor oj — —10°, and for a range of induced strains 
a. The curves are marked also by the voltages as we infer 
a simple linear relationship between a and the voltage to 
facilitate comparison with the experimental paper. 

The anisotropy behaviour shown in Fig. [221 can be de- 
scribed as a smooth rotation of the global energy mini- 
mum upon increase of a rather than the " scissors" effect 
shown in Fig. [5] in Sec. IIIIBI The total induced strain 
now contains both components e'' and e" as written in 
Eq. (pS)) . The uniaxial basis strain e" present due to the 
tilt of the stressor diminishes significantly one of the lo- 
cal minima typically occurring because of interplay of a 
positive cubic and a small uniaxial anisotropy component 
along a crystal diagonal. The remainder of the weaker 
local minimum is observed theoretically for a correspond- 
ing to voltages « —100 V when the component of the 
induced strain and the "intrinsic" shear strain compen- 
sate each other. One would expect domination of cu- 
bic anisotropy with two equivalent local minima close to 
the main crystal axes if the stressor had purely diago- 
nal alignment. The e" component of the total strain of 
the tilted stressor makes the local minimum closer to the 



[010] axis less pronounced (marked by arrow in Fig. 

For completeness, we discuss the free energy phe- 
nomenological formula used in Ref. [3l| to describe the 
in-plane angular dependence of the induced anisotropy. 
The decomposition of the total induced strain in Eq. ([26]) 
into the strain basis introduced in Eqs. ([MTTjl is not con- 
sidered in that work. Instead, the induced anisotropy is 
described by a single uniaxial term Kq sm'^{^ — fl) added 
to the phenomenological formula rather than terms with 
coefficients K[iio] and i^[iio] from Eq. ([55)) . Effectively, 
this corresponds to a change of variables from -/^[iio] 
and i^[iio] to Kq and fl. The angle O is measured 
from the [110] axis and it rotates the additional uniax- 
ial anisotropy term so that it describes the effect due to 
the tilted stressor. One may assume coUinearity of the 
resulting anisotropy component with the principal elon- 
gation direction of the piezo. However, this simple situ- 
ation is observed both theoretically and experimentally 
only when the stressor is aligned with the main crys- 
tal axes or diagonals. The missalignment for arbitrary 
orientation of the induced strain is due to the under- 
lying cubic symmetry of the system incorporated into 
our microscopic band structure calculation in the form 
of the band parameters 72, 73, 02, and 03. It has been 
explained in Sec. IIII CI that the coUinearity of the in- 
plane strain and corresponding anisotropy occurs only 
for the strains or e" (see Eqs. ([TUl) and ([TT]) ). For 
any other stressor orientation, ft ^ uj, which is reflected 
on the level of the anisotropy functions by the inequality, 
9[ioo] (2^7 p) 7^ ^[iio] i^iP)- It expresses the difference in the 
effect on magnetic anisotropy between straining the lat- 
tice along the main crystal axis and along the diagonals 
(see Eq. ([T7|) in Sec. HCTCl) . 

The transformation from variables K[iiq^{x,p,lu) and 
K[iio]{x,p,uj) to Kq_{x,p,uj) and Q,{x,p,uj) in the phe- 
nomenological formula in Eq. 
reads: 



for -7r/2 < uj < tt/2 



Fu{M) = if[no](c^)sin2i^ + i^[ioo](c^)sin2(i^-f 7r/4) 
= Knsm^[^-n), (40) 

where: 



1 



Vl(x,v,uj) — — arctan , 
^ ^ 2 \ K\ 



[100] 



110] 



(41) 



Kfi{x,p,uj) = -i^[iio] cos2f2 -I- ifjioo] sin20. 



Considering the approximate relation q[ioo] — 0.43g[iio] 
the formulae in Eq. ([4T|) simplify to: 



^ 1 , /'7[ioo](a;,p)sm2cj\ 

n{x,p,uj) = -arctan — — (42) 

2 \q[iio]{x,P)cos2LjJ 

— — arctan (0.43 tan 2uj) , 

q(i{x,p,uj) = q[iiQ]{x,p)cos2u;cos2fl + 

-|-0.43g[iio] {x,p) sin 2uj sin 2Q, 

where Kq = qQ_{a — a') /2. (The same transformation 
of variables can be used in case of strains induced along 
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arbitrary in-plane direction by relaxation in a narrow bar 
(see Eqs. P3II34|) ). Then we obtain = -gneop/2.) 

Note that in the representation of Fu{M) via i^[iio] 
and i^[iio] the dependence on lo can be simply factored 
out and the dependence on x and p is contained only in 
the functions functions (7[iio] and (/[loo] • For our general 
discussion presented in this paper it is therefore the more 
convenient form than Fu{M) expressed via Kn and fJ. 

We conclude that the in-plane alignment of the easy 
axis in patterned or piezo-stressed samples can be de- 
scribed on a semi-quantitative level by our modelling sim- 
ilarly to the bare (Ga,Mn)As epilayers. 



V. SUMMARY 

The objective of this work was to critically and thor- 
oughly inspect the efficiency of a widely used effective 
Hamiltonian model in predicting the magneto-crystalline 
anisotropics in (Ga,Mn)As. We have provided overview 
of the calculated anisotropies which show a rich phe- 
nomenology as a function of Mn concentration, hole den- 
sity, temperature and lattice strains, and compared it to a 
wide range of experimental works on the level of the mag- 
netic easy axis direction and on the level of anisotropy 
fields. The large amount of analysed results compensates 
for the common uncertainty in sample parameters as- 
sumed in experiment and allowed us to make systematic 
comparisons between theory and experiment on the level 
of trends as a function of various tunable parameters. 
Generically, we find this type of comparison between the- 
ory and experiment in diluted magnetic semiconductors 
much more meaningful than addressing isolated samples, 
given the complexity of these systems and inability of 
any theoretical approach applied to date to fully quan- 
titatively describe magnetism in these random-moment 
semiconducting ferromagnets. 

In Sec. in] we introduced the mean-field model used 
throughout the study, estimated the relative strength of 
the shape anisotropy, and discussed the correspondence 
of the shear strain, modelling the broken in-plane sym- 
metry measured in most (Ga,Mn)As epilayers, with a 
microscopic symmetry breaking mechanism. 

In Sec. mil we focused on modelling and experiments 
in bare unpatterned epilayers. The in-plane and out-of- 
plane magnetisation alignment was studied. For com- 
pressively strained samples the generally assumed in- 
plane anisotropy is found to be complemented by regions 
of out-of-plane anisotropy at low hole densities and low 
temperatures. This observation is corroborated by avail- 
able experimental data showing in-plane anisotropy in 
most of the studied epilayers but also the occurrence of 
the out-of-plane easy axis in materials with high hole 
compensation. At the same time, the model predicts 
out-of-plane easy axis for high hole densities at all Mn 
concentrations which has yet not been observed experi- 
mentally. 

Next, the competition of cubic and uniaxial in-plane 



anisotropy components was investigated. Wealth of ex- 
perimentally observed easy axis transitions driven by 
change of temperature or hole density finds correspond- 
ing simulated behaviour. The following general trend 
is observed in most samples: at low temperatures the 
easy axes are aligned close to the main crystal axes, 
while at high temperatures there is always diagonal align- 
ment. This trend is in good agreement with our calcu- 
lation, however, at low hole densities the calculated and 
measured easy axis transitions are more consistent than 
at higher hole densities where the measured phenomena 
match the predictions assuming hole densities typically 
a factor of two lower than in the experiment. 

We next introduced anisotropy fields corresponding 
to the crystal symmetry and to three distinct uniaxial 
strains. We extracted these anisotropy fields from the 
calculated data and found their dependence on mate- 
rial parameters. We observed linear dependence of the 
uniaxial anisotropy fields on the corresponding strains. 
Analysing experiments which determine the anisotropy 
fields from FMR, AMR or SQUID measurements allowed 
for detailed comparison of the cubic anisotropy compo- 
nent and two uniaxial anisotropy components (due to 
growth and the [110/[110] symmetry breaking). The 
measured and calculated anisotropy fields are of the same 
order of magnitude (~ 10^ — 10'^ Oe) in most samples. 

Finally, in Sec. IIVI we investigated structures where 
the post-growth patterning or piezo stressing was used 
to induce additional strains along any in-plane direction. 
The interplay of the intrinsic and induced anisotropies 
was studied. We discussed the procedure for obtaining 
the strain Hamiltonian from the parameters describing 
the experimental setup and a finite element solver was 
employed to find the inhomogeneous lattice relaxation in 
the patterned epilayers. Induced anisotropies were cal- 
culated directly using the total strain tensor. Alterna- 
tively, we also introduced a decomposition of the total 
strain matrix for any of the studied materials and device 
configurations into three basis strains and their additive 
effect on the total anisotropy. We found an overall semi- 
quantitative agreement of theory and experiment on the 
level of easy axis reorientations due to induced strains. 

The limitations of the theory approach employed in 
this paper have been thoroughly discussed in Ref. 0]. 
The model, which treats disorder in the virtual crystal 
approximation and magnetic interactions on the mean- 
field level is expected to be most reliable at lower tem- 
peratures and in the (Ga,Mn)As materials with metallic 
conductivity. We have shown that despite the limita- 
tions, the model captures on a semi-quantitative level 
most of the rich phenomenology of the magnetocrys- 
talline anisotropies observed in (Ga,Mn)As epilayers and 
microdevices over a wide parameter range. We hope 
that our work will provide a useful guidance for future 
studies of magnetic and magnetotransport phenomena in 
(Ga,Mn)As based systems in which magnetocrystalline 
anisotropies play an important role. 
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APPENDIX A: SYMMETRIES OF THE 
KOHN-LUTTINGER HAMILTONIAN 



The states at the top of the valence band have p-like 
character and can be represented by the l—l orbital mo- 
mentum eigenstates \l,mi). In the basis of combinations 
of orbital angular momentum eigenstates: 



l^> = ^(|1-1) + |11)), 



l^> - -^(|1-1)-|11)), 
\Z) = |10) 



(Al) 



Different representations of the six-band Kohn- 
Luttinger Hamiltonians are used in literature. Here, the 



the Kohn-Luttinger Hamiltonian for systems with no 
spin-orbit coupling can be written as: 



where 



hp 



e^, + Akl + B{kl + kl) 

Chy h'X 



A = 



B 



Ckx ky 

„2 I r>(u2 I U2\ 



ey+Ak:^ + B{ki + ki) 
Ckzky 



4 E 

TT7f ^ ^ 



\{x\pM' 



\{X\Py\^)\' 



2mo mi 



E 



Ckxkz 
Ckykz 
+ Akl + B{kl 



{X\px\i){t\py\Y) + {X\py\i){i\p,\Y) 



i^{X,Y,Z} 



ei 



(A2) 



(A3) 
(A4) 
(A5) 



and e„ is the energy of the valence band p-orbitals. 

The simple form is due to the symmetry of the zinc- 
blende crystal structure. The summation in elements 
A, B, C runs only through the Fi and r4 states of 
the conduction band as other levels are excluded by the 
matrix element theorem combined with the tetrahedron 
symmetry.^^ The only nonzero momentum operator ex- 
pectation values with neighbouring states are: 

{X\py\Ti{z)) = {Y\p,\Ti{x)) = {Z\p.,\Vi{y)) 

{X\p,\T,) = {Y\py\T,) = (Zb,|Fi). (A6) 

Due to the reflection symmetry with respect to the (110) 
planes it holds alsoi^ 

{X\py\T4z)) = {Y\p,\n{z))) (A7) 



If the tetrahedral symmetry of the GaAs lattice is bro- 
ken by potential V — xy^ as described in Sec. Ill Al the 
states Fi and ^iiz) of the conduction band considered 
in the summation in Eq. ([4]) are mixed, whereas states 
F4(a;) and F4(y) are left unchanged. In the perturbed ba- 
sis aTi+PT4{z), — /3T I + aT 4{z) , F4(a;), F4(j/) we obtain 
terms containing the parameter D in the Hamiltonian 
Hkp (See Eq. ^ in Sec. Ill AH A weak local potential V, 
a » P was assumed so terms of quadratic and higher 
order dependence on V could be neglected. Therefore the 
expression for parameters A, B, and C does not change. 
Using Eqs. (|A6P and (|A7P allows also for a compact ex- 
pression of the parameter D: 



D = C(Ab^|F4(z))(Fl|p,|A), (AS) 
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c 



1 



1 



-(ecl + A) e„-(ec4-A)_ 

where ed and are the energies of the conduction band 
Fi and r4 states, respectively. The small energy A is 
quadratically dependent on the size of the potential V 
but we include it to express the shift of the perturbed 
energy levels. 

To include spin-orbit coupling we use the basis formed 
by combinations of orbital angular momentum: 



|1> 
|2> 
|3> 
|4) 
|5) 



\J 
\J 
\J 
\J 



3/2, 771 j 

3/2, 
3/2, m^- 

3/2, 77lj 



3/2) 
-1/2) 
1/2) 
-3/2) 



(A9) 



|j = l/2,m, =.1/2) 
|6) s |j = l/2,m, = -l/2) 

The spin-orbit correction to the 6-band Hamiltonian is 
diagonal in this basis and can be parametrised only by 
a single parameter Aso-^ The 6-band Kohn-Luttinger 
Hamiltonian in the representation of vectors (|A9p reads: 
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bV3 
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(AlO) 

The 4-band Hamiltonian is highlighted. The Kohn- 
Luttinger eigen-energies are hole energies (measured 
down from the top of the valence band). The matrix 
elements of Hkl are listed in Ref. [5| . Here we focus on 
the modification of these elements due to incorporating 
the microscopic potential V = xy^: 

'Hhh 
b 



h 

^ [(71 +I2)ikl + kl) + (71 - 2j2)kl + G^ik^ky 



I2){kl + kl) + (71 -I- 272)fc2 + 2j^k,ky] 



2m 



^1 



^kl) 
iky) 



^^Akxky] 



V3?i 
2m 
_V2^ 
2m 



liikl 
lzk^(h 

b2{kl 

[j2i2kl - kl - kl) ~ 2j4k.ky 



2i{"fikxky + '^{kl 



(All) 



where we neglect the higher order effect of broken sym- 
metry on standard Luttinger parameters: 

2too 



71 

72 

73 



^2 {A + 2B) 



3h^ 



(A 
C 



B) 



and add a new parameter: 



74 



2mo 



D. 



(A12) 



APPENDIX B: LATTICE STRAINS AND 
MICROSCOPIC POTENTIAL 



We incorporate the lattice strain into the k • p the- 
ory following Ref. [25*1, which shows that Hkl and the 
6-band strain Hamiltonian have the same structure given 
in Eq. (jA10[) . Components of the strain tensor ea/3 in- 
troduced in Eq. ([3]) play role of the fc- vector components. 
The replacements in of matrix elements of Eqs. (|Alip (or 
rather of Eq. (A9) in Ref. 0) read: 



kakh 



Sq/3 



(Bl) 



2too 



71 



ai, 



7,mn 1^ 2 ' 



2mo 



73 



as 

2V3' 



where ai, 02, and 03 are the elastic constants. Their val- 
ues are different to Luttinger parameters 71, 72, and 73 as 
they originate from the first order momentum operator 
perturbation due to strain and second order perturba- 
tion treatment of the k • p term, respectively. The strain 
Hamiltonian has the following elements (in the hole pic- 
ture): 



'^hh — 


- (oi + y ) [exx + eyy) - 


(ai - a2)e 


njS 

n-ih = 


— ^fli {Gxx + ^yy) ^ 


(ai + a2)e 


Tilo = 


ai{Gxx ^yy ^zz) 






-03(6^2; - ie^_y) 






~^V^(Cyy ^xx) ~^ ia^Cxy 




d" = 


V2 

—a2 {2e.z - {exx + Cyy)) 





(B2) 



Now we compare the effect of microscopic symmetry 
breaking described by including the 74 dependent terms 
into the Hamiltonian Hkl to the effect of a uniform lat- 
tice strain incorporated as Tistr with matrix elements 
given in Eqs. (|B2[) . First, we write the strain Hamil- 
tonian Tistr as a sum of a contribution corresponding 
to the in-plane shear strain along the [110] axis and the 
growth strain introduced in Sec. IHI CI by Eq. ^ and 
Eq. pU)) . respectively. Their magnitudes are denoted by 
Cxy and 

^xx — ^yy 

= cq. Then we write the correc- 
tion Hv — T^KL — T^KL to the 6-band Kohn-Luttinger 
Hamiltonian due to the microscopic potential V = xy^ 
breaking the tetrahedral symmetry of the crystal as a 
sum of terms with different dependence on the in-plane 
direction of the fc- vector: 
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(B3) 



,(B4) 



where cn, Ci2 are the elastic moduh. ^'^i^^ 

Note that by resetting the reference energy in Eq. (|B4|) 

by subtracting S-^j^k^ky from the Hamiltonian Hv and 
the following substitutions: 

V^h^ ,,2 , ,2\ 

Cll + 2ci2 Ti^ , . 

-0260 -f4kxky (B5) 

Cll mo 

we can identify the two components of the strain Hamil- 
tonian in Eq. (|B3p with the two components of the Hamil- 
tonian Hv in Eq. (jB4p . The important difference, how- 
ever, is the dependence on fc-vector in case of Tiy- The 
first term of Tiv depends on the magnitude of the fc- 
vector, not on its in-plane orientation. The second term 
of Hv has the same structure as the second term of Hstr 
(which incorporates the effect of the growth strain) , how- 
ever, it does depend on the in-plane direction of the fc- 
vector so it generates a uniaxial in-plane anisotropy com- 
ponent that contributes to the energy profile (shown in 
Fig. [1]) similarly to the first term of Eq. (|B4p (contrary to 
the negligible uniaxial in-plane anisotropics correspond- 
ing to the growth strain). 

APPENDIX C: CUBIC ANISOTROPY TERMS 

The angular dependence of the magnetocrystalline 
anisotropy energy can be approximated by a series of 
terms of distinct symmetry. In Sec. lIHCl we introduced 
a simple phenomenological formula consisting of the low 
order terms of the cubic and uniaxial symmetry. Here we 
explain the choice of the independent cubic terms. 

We write the terms using the components of the 
magnetisation unit vector M: Ux = cos ^ sin 0, Uy — 
sin(/)sin0, = cos 9, where our angles 9 and (j) are mea- 
sured from the [001] and [100] axis, respectively. The cu- 
bic symmetry requires invariance under permutation of 



the coordinate indices x, y, and z. The simplest term sat- 
isfying the condition is equal to unity: n^+riy + nl = 1. 
The first order cubic term can be derived from its second 
power: 

{nl + nl + nlf^ (CI) 

n r 2 2 , 22, 2 2\ 

= 2 [n^Uy + n^n^ + Uyin^) 

We obtained two lowest order cubic terms which are mu- 
tually dependent. Therefore it is enough to choose only 
one of them. In case of Eq. the lowest order cu- 
bic anisotropy term reads: Kd {n^n^ + "^"z + ^^^y)^ 
where Kd is an energy coefficient. 

The second order term can be derived from the first 
order term: 

{nlnl + nlnl + nlnl) {nl + 71^ + nj) = (C2) 

42, 42i 24, 42, 24, 24, 

= nxny + n^n^ + n^Uy + UyU^ + n^n^ + UyU^ + 

I 2 2 2 

+ n^UyU^. 

The two second order terms and the first order term are 
mutually dependent. Again, only one term describes fully 
the second order component of the cubic anisotropy. We 
choose Kc2 {n-xiT-y^l) to be included into our approximate 
formula in Eq. ([fe)) . 

The independent third order term is derived as follows: 

{nlnl + nlnl + n^n^) {nl + + nlf ^ (C3) 

o/'422i 242, 2 2 A\ 

= ^(nxnyn.+n^nyn^ + n^nynj 
+ "^intn^ + nlnt + n^ni) 

+ nlnl + nlnl + nln^ + nlnl + nlnl + n^n^. 

Note that the first part of the product is proportional to 
the second order cubic term. Again, we can choose one 
of the two dependent third order terms. This derivation 
procedure can be continued but fitting our microscopic 
data to the phenomenological formula yields a negligible 
magnitude even for the third order term coefficients. 
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APPENDIX D: USED CONSTANTS 



Let us list all the material parameters used in our codes 
for (Ga,Mn)As: 



71 


72 


73 


6.85 


2.1 


2.9 


ai[eV] 


a2[oV] 


a3[cV] 


-1.16 


-2.0 


-4.8 


cii[GPa] 


ci2[GPa] 


a;c[nm] 


12.21 


5.66 


0.565 




Jpd [eVnm^] 




0.341 


0.055 
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